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Executive  Summary 


Digital  cameras  use  a  focal  plane  array  (FPA)  to  convert  light  energy  into  electrical  energy 
that  can  be  processed  and  stored  in  the  memory  chip.  The  FPA  is  composed  of  pixels  or  chips  of 
electro-optical  material  on  a  planar  surface.  Placing  the  sensors  (pixels)  on  a  plane  surface  can 
lead  to  off-axis  (spherical)  aberrations  and  can  limit  the  field  of  view  (FOV)  of  optical  systems 
unless  exotic  optics  are  used.  To  correct  for  these  distortions,  designers  have  to  use  additional 
optical  elements,  which  complicate  the  design  of  the  optics  and  increase  the  cost  of  the  cameras. 

Curved  FPAs  with  simple  spherical  lenses  would  provide  improved  performance  (vs.  the 
commonly  used  planar  FPA).  These  curved  FPAs  would  provide  a  large  FOV  with  better  resolu¬ 
tion  off  axis,  would  require  fewer  lenses,  and  would  eliminate  the  need  for  image  post 
processing.  The  systems  that  employ  these  curved  FPAs  have  superior  optical  properties,  but  the 
curved  FPAs  are  more  difficult  to  manufacture.  To  meet  this  challenge,  the  Defense  Advanced 
Research  Projects  Agency/Microsystems  Technology  Office  (DARPA/MTO)  has  instituted  a 
program  called  the  Hemispheric  Array  Detector  for  Imaging  (HARDI).  The  program’s  goal  is  to 
develop  curved  FPA  technology  and  combine  it  with  appropriate  lens  systems  to  enhance  mili¬ 
tary  capability.  To  accomplish  this  goal  of  creating  a  curved  FPA,  the  HARDI  program  is 
exploiting  the  properties  of  organic  and  hybrid  organic/inorganic  semiconductor  materials. 

In  support  of  the  HARDI  program,  the  Institute  for  Defense  Analyses  (IDA)  performed 
technical  assessments  and  provided  planning  assistance.  IDA  has  identified  promising  applica¬ 
tions  of  spherical  FPAs  that  could  potentially  be  exploited  for  defense-related  applications. 
Curved  FPAs  remove  the  need  for  complex  lens  systems  to  correct  off-axis  distortion,  and  this 
improvement  alone  reduces  the  weight,  size,  and  complexity  of  lens  systems.  From  the  multiple 
specific  applications  of  curved  FPAs  identified  by  IDA  (cameras  mounted  on  small  robots, 
miniature  unmanned  aerial  vehicles  (UAVs),  and  small  surveillance  cameras),  HARDI  man¬ 
agement  decided  to  focus  first  on  the  cameras  used  on  small  robots  for  the  Advanced  Mine 
Detection  System  (AMDS)  Program.  HARDI  management  then  directed  IDA  to  provide  an  inde¬ 
pendent  analyses  of  these  proposed  systems,  including  appropriate  selection  of  lenses  and 
tradeoff  studies  of  pixel  size  and  count  vs.  camera  characteristics  to  obtain  optimum  conditions 
for  various  robot  camera  applications.  For  this  task,  IDA  developed  a  ray  tracing  code,  compared 
this  code  to  published  results  and  simple  analytical  closed-form  solutions,  and  began  using  it  to 
study  applications  of  interest  to  HARDI. 
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1.  Introduction 


A.  Background 

Digital  cameras  use  a  focal  plane  array  (FPA)  to  convert  light  energy  into  electrical 
energy  that  can  be  processed  and  stored.  The  FPA  is  a  two-dimensional  (2-D)  array  of 
photodetectors  (or  pixels)  fabricated  on  an  electro-optical  material.  Used  in  this  way,  the 
term  pixel  refers  to  the  actual  physical  detector  in  a  camera.  In  digital  imaging,  however, 
a  pixel  is  the  smallest  addressable  picture  element  in  raster  graphics — the  smallest  unit  of 
picture  that  can  be  controlled.  Figure  1-1  presents  both  concepts  of  a  pixel.  In  this  report, 
we  will  use  the  concept  of  pixel  represented  in  Figure  l-l(b). 


(a)  (b) 

Figure  1-1.  Two  Usages  of  the  Term  Pixel:  (a)  Pixel  Pattern 
as  Used  in  Imaging  and  (b)  Physical  Pixels  (Bottom)  Forming 
an  FPA,  With  Micro  Lenses  (Top)  and  Color  Filters  (Center) 

The  size  and  the  number  of  the  individual  pixels  affect  the  resolution  and  the  read¬ 
out  speed  of  the  camera.  The  more  pixels  a  digital  camera  has,  the  more  detail  it  can 
record  when  a  photo  is  taken.  Smaller  pixels  improve  resolution,  if  the  camera  resolution1 
is  not  optics  limited,  but  may  reduce  sensitivity  (i.e.,  need  more  time  for  required  energy 
deposition).  Modem  digital  cameras  contain  FPAs  that  have  pixel  counts  on  the  order  of 
megapixels.  For  example,  cameras  with  2  megapixels  are  becoming  obsolete,  cameras 
with  5  megapixels  are  in  decline  but  still  a  good  value,  and  cameras  with  10  megapixels 
are  in  the  mainstream. 


Resolution  of  systems  is  generally  described  by  the  modulation  transfer  function  (MTF),  which  is 
defined  and  discussed  at  the  beginning  of  Section  4. 
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In  current  digital  cameras,  a  design  constraint  of  optical  systems  is  that  the  image 
surface  (or  Petzval  surface)  must  be  planar  so  that  the  image  can  be  recorded  on  a  planar 
silicon  FPA.  This  requirement  leads  to  off-axis  aberrations  including  spherical,  astigma¬ 
tism,  field  curvature,  and  coma.  To  correct  for  these  distortions,  designers  use  additional 
optical  elements,  which  complicate  the  design  of  the  optics  and  increase  the  cost  of  the 
cameras.  Figure  1-2  shows  a  modem  camera  with  the  complex  lens  system  and  the  flat 
FPA. 


Lens  System  Focal  Plane  Array 


Figure  1-2.  Cutout  of  a  Camera  Showing  the  Lens  System  and  the  Position  of  the  FPA 

To  avoid  the  complexity  of  lens  systems  shown  in  Figure  1-2,  the  use  of  curved 
FPAs  with  simple  spherical  lenses  has  been  proposed  and  is  being  developed  for  military 
applications  by  the  Defense  Advanced  Research  Projects  Agency/Microsystems  Tech¬ 
nology  Office  (DARPA/MTO)  as  part  of  the  Hemispheric  Array  Detector  for  Imaging 
(HARDI)  program.  Figure  1-3  shows  ray  tracing  results  for  three  optical/FPA  systems 
and  illustrates  the  benefits  of  curved  FPAs  over  the  flat-plane  FPAs. 


(a)  (b)  (c) 

Figure  1-3.  Comparison  of  Ray  Tracings  for  Different  Systems 
Note  for  Figure  1-3:  (a)  Single  element  lens  with  a  planar  FPA,  (b)  Ball  Lens  with  a  curved  FPA,  and  (c) 
Three-element  lens  system  with  a  planar  FPA. 


Joseph  Petzval  (1807-1891)  was  a  Hungarian  mathematician,  inventor,  and  physicist.  He  is  best  known 
for  his  work  in  optics.  Petzval  is  considered  one  of  the  main  founders  of  geometrical  optics,  modem 
photography,  and  cinematography.  Among  his  inventions  are  the  Petzval  portrait  lens  and  opera  glasses, 
both  still  in  common  use  today. 


1-2 


Ray  tracing  results  for  a  single-element  plano-convex  lens  (a  commonly  used  simple 
lens)  on  a  planar  FPA  indicate  some  obvious  problems  (see  Figure  l-3(a)).  First,  notice 
that  the  off-axis  rays  do  not  focus  on  the  FPA  beyond  the  center  of  the  FPA.  This  lack  of 
focus  causes  blurring  of  point  spots  and  reduces  the  resolution  of  the  system.  Reduction 
of  relative  illumination,  referred  to  as  fall-off,  at  the  edge  of  pictures  also  occurs  because 
the  focal  length  of  the  ray  does  not  match  the  location  of  the  planar  FPA.  This  issue  of 
ray  spreading  affects  proper  data  collection  and  will  be  discussed  in  more  detail  in  Chap¬ 
ter  4  when  we  define  the  point  spread  functions  (PSFs)  and  modulation  transfer  functions 
(MTFs)  and  use  them  to  compare  different  lens  systems. 

The  system  with  the  spherical  FPA  and  ball  lens  with  the  same  center  point 
improves  ray  focusing  and  light  illumination  at  different  angles.  All  the  rays  are  axial  and 
pass  close  to  the  center  of  the  lens  and  the  FPA  sphere  because  they  are  constrained  by 
the  aperture  in  the  center  of  the  sphere  (see  Figure  l-3(b)).  The  symmetry  of  the  system 
when  the  rays  are  constrained  to  pass  through  the  center  of  the  spherical  lens  and  the 
spherical  FPA  allows  for  a  wide  field-of-view  (FOY). 

To  maintain  good  image  quality  while  using  a  flat  FPA,  multiple  lens  systems  have 
to  be  designed  to  correct  distortions.  Figure  l-3(c)  illustrates  the  ray  tracing  results  for  a 
lens  system  using  three  lenses  and  terminating  on  a  flat  FPA.  Although  this  shows  an 
improvement  over  Figure  1-3  (a),  increasing  the  FOV  can  still  be  a  problem. 

B.  HARDI  Program/Institute  for  Defense  Analyses  (IDA)  Task 

1.  HARDI  Program 

State-of-the-art  cameras  could  be  improved  by  increasing  their  FOVs  with  uniform 
illumination  over  all  angles  and  decreasing  their  weight  and  cost.  A  simple  curved  FPA 
would  provide  this  improvement. 

Curved  FPAs  have  been  suggested  as  having  superior  optical  properties,  but  they 
are  more  difficult  to  manufacture  and  have  been  restricted  to  large  radii  of  curvature  until 
recently.  New  technologies  have  allowed  the  manufacture  of  curved  FPAs  with  diameters 
as  small  as  1  cm. 

To  develop  curved  FPA  technology  and  combine  it  with  appropriate  lens  systems  to 
improve  military  capability,  DARPA/MTO  has  put  in  place  the  HARDI  program.  The 
objective  of  this  program  is  to  exploit  properties  of  organic  and  hybrid  organic/inorganic 
semiconductor  materials  to  create  a  curved  FPA.  This  curved  FPA  will  provide  a  large 
FOV  that  has  better  resolution  off  axis,  will  require  fewer  lenses,  and  will  eliminate  the 
need  for  image  post  processing. 

The  first  two  phases  of  the  HARDI  program  are  approaching  completion.  The  intent 
of  the  third  phase  is  to  provide  a  useful  device  for  an  important  military  application.  In 
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other  words,  DARPA  needs  to  focus  on  a  real  military  problem  and  optimize  the  HARDI 
system  to  meet  the  needs  of  the  military.  To  do  this,  we  need  to  know  the  current  capabil¬ 
ities  of  military  systems  that  might  be  improved. 

2.  IDA  Task 

IDA  performed  technical  assessments  and  provided  planning  assistance  to  the 
HARDI  program  by  evaluating  the  three  different  approaches  to  implementing  curved 
FPAs  into  military  systems.  The  task  statement  included  these  areas  for  investigation: 

•  Identify  military  systems  that  can  be  simplified  and  substantially  improved  by 
combining  curved  FPAs  with  simple,  small  lightweight  commercial  lenses; 

•  Identify  potential  innovative  applications  by  combining  curved  FPAs  with  sim¬ 
ple  commercial  lenses  (e.g.,  determining  the  feasibility  of  using  rail-mounted, 
moving  curved  FPAs  to  provide  passive  three-dimensional  (3-D)  target  imaging: 
detection,  location,  and  identification);  and 

•  Identify  a  potential  application  for  using  Gradient  Index  of  Refraction  (GRIN) 
lenses  with  curved  FPAs  to  provide  wide  FOV  optical  systems  (e.g.,  gun  and 
weapon  sights). 

In  support  of  the  HARDI  program,  IDA  has  initially  focused  on  the  first  area  and 
has  identified  promising  applications  of  spherical  FPAs  that  could  potentially  be 
exploited  for  defense-related  applications.  This  discussion  is  included  in  Chapter  2  of  this 
report.  In  the  past,  flat  image  planes  increased  the  mass  and  complexity  of  camera  optics. 
Curved  FPAs — like  the  eye’s  retina — alleviate  this  problem  by  removing  the  need  for 
complex  lens  systems  to  correct  off-axis  distortion.  This  improvement  alone  reduces  the 
weight,  size,  and  complexity  of  lens  systems  in  military  applications. 

From  the  multiple  specific  applications  of  curved  FPAs  identified  by  IDA,  HARDI 
management  has  decided  to  focus  first  on  cameras  used  on  small  robots  and  directed  IDA 
to  provide  information  about  these  cameras  (e.g.,  their  technical  specifications  and  oper¬ 
ational  capabilities).  These  robots  are  used  in  the  Advanced  Mine  Detection  System 
(AMDS)  and  Explosive  Ordnance  Disposal  (EOD)  programs.  IDA  obtained  technical 
information  about  some  of  these  cameras  and  investigated  their  operational  capabilities  in 
the  field.  This  information  is  provided  in  Chapter  3  of  this  report. 

IDA  was  then  asked  to  provide  independent  analyses  of  these  proposed  systems. 
These  analyses  included  appropriate  selection  of  lenses  and  tradeoff  studies  of  pixel  size 
and  count  vs.  camera  characteristics  such  as  resolution,  FOV,  and  so  forth  to  obtain  opti¬ 
mum  conditions  for  various  robot  camera  applications.  For  this  task,  IDA  developed  a  ray 
tracing  code  (see  Appendix  A  of  this  document),  compared  the  code  to  published  results 
and  simple  analytical  closed-form  solutions,  and  began  using  it  to  study  systems  of  inter¬ 
est  to  HARDI.  This  activity  is  discussed  in  Chapter  4  of  this  report. 
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2.  Potential  Military 
Applications  of  Curved  FPAs 

Chapter  2  presents  the  HARDI  Follow-Up  Meeting  briefing  that  IDA  presented  to 
DARPA/MTO  on  July  21,2010. 
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HARDI  Follow-Up  Meeting 


O 


ru 

CD 


O 

CD 


IDA 


4-  £ 

O  E 

co  CD 


ru  — 

• 

£  E 
o 

o  c 

Q_  — 
CD 

CD  LL. 

- - 

CD  CO 

3  > 

IS 

LU  (T3 


V 

> 

O  co 
s-  CD 

a  CO 

.§  S 

>»  +-» 

"ro  cud 

CD 


"O  — 
CD  Q 


</>  — 

c  I 

CD  co 

■O  <u 

CD 

:—  Q_ 

Si 

E  ” 

<o  +f 
o  5 

c  ^ 
CD  Q- 

U  LL. 


CO  Z5  +J 

£  U  CD 

a>  bo  g 

+-J  C  g 

co  •—  C 

>  c  c 

CO  •—  .— 

>  “2  > 

M—  F  M— 

+-»  o  *-* 

C  u  C 

CD  ^  CD 

2  -Q  2 


g  C 

co  3 

<  ° 

£| 

■g  e 

£  cuo 

u 

cuo  3 

C  4— 

c  ° 

s  ^ 

i! 

u  00 

>  CD 

-Q  £ 

CO 

§  a? 

•43  CD 

CD  ^ 
CJ  co 
—  CD 
Q_  co 

ro  .1 

>  ™ 

5  <u 

o  E 


c  “■ 

■S  -i 

—  co 


co  CD  "o 

^  ^  " 's 

CD  ^  Cg 
£  <=  —  2 

•-  o  +;  q 

“O  ’J3  c 

CD  S  £  2 

CD  . —  ^  in 


■M  c 
CD  CD 
>  “D 


S..2 

CD  "fo 
“D  U 

>  O 


O 

o  +3 

■m  u 

co  CD 

S:  3L 

“O  Cuo 

CD  c 

£  'cuo 

D  (D 

u  E 

cuo  — 

C  -t— ’ 
•-  CD 

>  cuo 

O  cT 
C  (D 


o  g  w 

X  CD  -£ 
CD  "O  CUO 

II  1 

c  “■  g. 

(D  O  fD 
“  ^  (D 
^  co  ^ 

g*n  ™ 

._  <D  c 
co  >  15 
3  1-  cuo 
s-  3  . 

£  £  b? 

.2  3  „ 

IS  »i  £ 

U  O)  CD 

1=  V)  +-» 

Q.  C  ^ 

S-®  5r 

(D  _ 

TZ  “a*"  fD 

(D  ^  (j 
+i  ^  -3 
^  c  n 

£  s.  o 

2  c  >■ 

£.2  6 

.>*  +■»  LL, 
M—  M  — - 

CD  T3  (D 

2  0£  > 


^j- 

<N 


(/) 

<U 

> 


u 


QJ 


<U 

Q) 


c  >• 

CD  CuO 

E  o 

cu  o 
>  c 
o  _c 

Q-  oj 

E  I- 

£  CU 

to  CD 
CD  U 

'§>.§ 

—  u 

o  on 
c  ^ 

(—  to 
u  CD 
CD  £ 

+-»  ji 

—  CU 

.E  c  ^ 

■p  <  tt 

CD  CD  ^3 
■m  co  co 
O  C  _ 

Q- <y  q 


■+T  CD 

°  Q  c/) 


CU  Q) 
U  +-» 


O 
co 

«  c  :> 
^  CZ-Q 


CO 


CO 

(D 

'cuO 

O 


CO 


CD  -3 
CD  — 
■M  O 
“O  co 
CD 

+j  * 

aJ  Q 

■53  2= 

co  <£ 


i  £ 

CO  LJ_ 

E  T5 

aj  cu 

s  £ 

o  3 

5-  U 
Q_  _ 
■  .  CU 


CD 


CD 


=5  O 
<->  Q. 


■?od  5 


O  o 

C  ^ 
O  CD 
h-  co 
CU  ^ 
.y  CU 


CD 


CU 

"D 


co 

&_ 

CD 

~o 

CU 

CD 


CJ 

CD 

B 

o 

i_ 

Q_ 
■  • 
co 


CO 

CO 

3 

u 

CO 


CO 

-M 

u 

cu 

-M 

C  10 

o  t 

CD  CD 
Q_ 
X 
CD 

■v 

to 

o 

■M 
CJ 
CU 


o 

c 

o 

B 

■M 

CU 

CJ 


2  ra 


a>  5 

2  8 


CO 

+-» 

U 

cu 


o 

u 


CO 

CD 

u 

cu 

o 

L_ 

Q_ 

Q_ 

CU 

CUO 


£  +3 

O  Q_ 


CO 

co 

3 

u 

CO 


CO 

tuO 


CD 

CD 


O 

CJ 


O 

B 

■M 

(U 


CU 

> 


in 

<N 


HARDI  =  Hemispheric  Array  Detector  for  Imaging 


Small  Robots  Used  in  the 
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Solution.  This  problem  could  be  attacked/helped  with  improved  off- 
axis  resolution 
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Miniature  Unmanned  Aerial  Vehicles  (UAVs) 
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I  he  Night  vision  ana  biectromc  sensors  Directorate  (NVbbU)  Air 
Systems  Division  (ASD)  collaborates  with  the  Army  Research 
Laboratory  (ARL)  on  imaging  systems  for  UAVs.  Need  to  establish 
a  contact 
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Check  to  see  how  wide  FOV  systems  with  flat  FPAs  compare  with  curved  FPA  systems 

Check  how  well  cell  phone  cameras  perform.  Compare  with  the  military  requirements.  These 
cameras  are  small  and  lightweight  and  have  low-cost  optics 


Objectives  of  task  presented 
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3.  Cameras  Used  With  Small  Robots 


Chapter  3  presents  the  HARDI  Follow-Up  Meeting  briefing  that  IDA  presented  to 
DARPA/MTO  on  October  5,  2010. 
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Erik  Rosen/Institute  for  Defense  Analyses  (IDA)  suggested  that  we  also  look  into  the  cameras  used  by  robots 
in  the  Explosive  Ordnance  Disposal  (EOD)  program,  which  faces  requirements  similar  to  those  in  the  AMDS 
program 


Production  of  Miniature  Robots  Used  by  AMDS 
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QinetiQ  North  America  (Talon) 
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size  used  on  another  similar  camera:  2.85  |jm 


Zoom  Camera  on  Talon  (1) 
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Zoom  Camera  on  Talon  (2) 


Tracks 


Christopher  Geyer,  (Lead  Scientist  -  Imaging 
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Dimensions:  50  x  57.5  x  81.8  mm  (W  x  H  x  D) 
Weight:  230  g 


Micro  Lens  -  Sensor  Structure 
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I  attended  QinetiQ  robot  demonstration  at  Quantico  Marine 
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Backup  Slides 
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Sony  Camera 


Sony  FCB-EX780B 

25x  Super  HAD  NTSC  Color  Block  Camera  with  External  Syn 

fvieir  product  brochure  in  ADOBE  PDF) 


This  new  FCB-EX780B  camera  is  an  evolution  in  security  dome  remote 
presence  and  traffic  monitoring  applications.  The  FCB-EX7S0B  is 
equipped  with  new  and  unique  surveillance  features  compared  to 
previous  FCB  models  such  as  an  E-flip  function  that  electronically  flips 
the  picture  for  correct  image  display  and  an  alarm  function  that  enables 
changes  to  be  detected  within  any  given  area  of  the  picture  In  addition, 
this  camera  features  an  improved  Privacy  Zone  Masking  function 
compared  to  previous  FCB  models  for  sophisticated  masking  privacy 
control  -  a  necessity  in  many  surveillance  applications.  Combining 
superb  picture  quality  that  you  expect  from  Sony  FCB  cameras  and  a 
variety  of  unique  and  convenient  features,  this  new  FCB-EX78QB 
camera  is  the  perfect  match  for  demanding  indoor  and  outdoor 
surveillance  applications. 

Adding  to  its  superb  flexibility  and  easy  operation  the  FCB-EX7S0B 
camera  incorporates  familiar  and  convenient  features  such  as  Spot  AE 
Auto  ICR  (IR  Cut  filter  removal)  quick  camera  control  via  a  high-speed 
serial  interface  (max  33.4  Kb/s),  and  various  customizable  settings. 

Features 

*  E:Flip  Function  (electronic  flip) 

*  25x  Optical  Zoom  /  12x  Digital  Zoom 

*  Image  stabilizer 

■  Auto  ICR  (IR  Cut  Filter  Removal) 

*  Alarm  Function 

*  Day/Night  Mode 

*  Picture  Freeze  during  zoom  focus,  preset  and  lens  initializing 

*  Key  switch  connector  (CN701)  and  DC/Video  Connector  (CN903) 

*  Privacy  Zone  Masking 

*  Minimum  Illumination  of  2.5  lux 

*  AE  Spot 

*  High-speed  serial  interface  (max.  33.4  kb/s) 

*  TTL  signal-level  control  (VI SCA  control) 

*  Various  factory  presets 

*  Intern  a  I/Extern  a  I  sync 

*  Lead -free  solder,  halogen  free  mounting  boards,  and  low  power  consumption 
(min.  1.6W  with  inactive  motors) 
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Sony  Camera  Specs 


Specifications 


Image  device 

1/6  type  Super  FIAD  CCD 

Number  of  Pixels 

Approximately  680.000  pixels 

Lens 

25Xzoom.  1=2.4  mm  (wide)  to  60  mm  (tele).  FI. 6 
to  F2.7 

Digital  Zoom 

12X  (30QXwith  optical  zoom) 

Angle  of  View 

45  degrees  (wide  end)  to  2.0  degrees  (tele  end) 

Min.  Object  Distance 

35  mm  (wide  end)  to  800  mm  (tele  end) 

Sync.  System 

Intemal/External  (V-Lock) 

Min.  Illumination 

2.5  lx  (typical)  (50  IRE) 

S/N  Ratio 

49  dB 

Electronic  Shutter 

1/1  to  1/10,000  s.  22  steps 

White  Balance 

Auto,  ATW,  Indoor.  Outdoor,  One-push,  Manual 

Gain 

Auto/Manual  (-3  to  28  dB.  2  dB  steps) 

AE  Control 

Auto:  Manual  Priority  mode,  Bright,  EV 
compensation,  Back-light  compensation 

EV  Compensation 

-10.5  to  +10.5  dB  (1.5  dB  steps) 

Back  Light  Compensation 

On/Off 

Privacy  Zone  Masking 

On/Off  (24  positions) 

Flicker  Cancel 

Auto 

Focusing  System 

Auto  (Sensitivity:  normal,  low).  One-push  AF, 
Manual,  Infinity,  Interval  AF.  Zoom  Trigger  AF 

Picture/Digital  Effects 

E-Flip,  Neg.  Art.  Black  &  White,  Mirror  Image 

Camera  Operation  Switch 

Zoom  tele,  Zoom  wide 

Video  Output: 

VBS:  1.0  Vp-p  (sync  negative),  Y/C  Output 

Camera  Control  Interface 

VISCA  fTTL  signal  level),  baud  rate:  9.6  Kb/s. 

19.2  Kb/s,  38.4  Kb/s,  Stop  bit:  1/2  selectable 

Storage  Temp 

-20’Cto  60aC  (-4"F  to  140aF) 

Operating  Temp 

Q’Cto  50aC  (32aF  to  122aF) 

Power  Consumption 

6  V  to  12  V  DC,  1.6  W  (motors  inactive J/2.5  W 
(motors  active) 

Weight:  approx. 

230  g  (8.1  oz) 

Dimensions 

50  x  57.5  x  81.8  mm  (2  x  2-3/8  x  3-1/4  inches) 
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Zoom  Camera  (QinetiQ) 
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Report  on  Micro  Lens  Sensor  Geometry 


1 /6-inch  380K/440K  Effective  Pixel 

Ultra  small  High -picture  Quality  Camera  CCDs 

ICX238AKE  (NTSC) 
ICX239AKE  (PAL) 

There  are  now  increasingly  strong  demands  for  even  further  miniatur¬ 
ization  and  improved  picture  quality  in  CCD  devices*  which  are  widely 
used  for  image  input. 

Sony  has  now  developed  two  new  high  picture  quality  ultrasmalt  CCD 
image  sensors  to  respond  to  these  demands. 

By  developing  the  industry's  smallest  unit  pixel,  a  pixel  with  a  horizon¬ 
tal  pixel  pitch  of  only  3.2  pm,  Sony  has  achieved  effective  pixel  counts 
of  380K  and  440K  in  the  ICX238AKE  and  ICX239AKE,  which  are  1/6- 
inch  optical  system  devices,  and  Sony  has  also  achieved  even  further 
reductions  in  device  power  consumption  in  these  products. 

These  features  can  contribute  to  increased  picture  quality,  further  min¬ 
iaturization,  and  reduced  power  consumption  in  digital  cameras. 


The  1CX238AKE  (NTSC)  and 
[CX239AKE  (PAL)  are  1/6- inch  opti¬ 
cal  system  color  CCD  image  sensors 
with  380K  and  440K  effective  pixels, 
respectively.  Despite  being  ultras  mall 
devices,  they  provide  the  high  level 
characteristics  required  by  camcorders, 
surveillance  cameras,  and  similar  prod¬ 
ucts.  Furthermore,  since  they  achieve  the 
low  power  consumption  and  other 
characteristics  required  by  notebook 
personal  computers,  PDAs,  and  other 
portable  information  terminals,  they  are 
optimal  for  use  in  a  wide  range  of 
applications, 

VOICE 


Despite  being  uttrasmall  devices, 
these  CCDs  achieve  excellent  pic¬ 
ture  quality  and  high  perfor¬ 
mance.  Thus  we  think  that  these 
are  optimal  CCDs  not  only  for 
current  products  that  use  CCDs, 
but  for  products  in  new  applica¬ 
tion  areas  that  desire  image  input 
functionality.  I  recommend  that 
you  will  look  into  using  the 
ICX23QAKE  or  ICX239AKE  In 
unique  new  products  that  can 
take  advantage  of  the  miniature 
size  and  low  power  consumption 
that  are  the  features  of  these 
devices. 


D  Further  Miniaturization  and 
Higher  Picture  Quality 

Due  to  the  development  of  fine  fabrica¬ 
tion  technology  that  can  create 
devices  with  a  horizontal  pixel  pitch  of 
3,2  ^fin,  which  corresponds  to  the 
industry's  smallest  unit  pixel,  Sony 
achieved  effective  pixel  counts  of  380 K 
and  44QK*  which  are  the  largest  pixels 
for  a  1 /6-inch  optical  system  device. 
Thus  these  devices  achieve  both  further 
miniaturization  and  higher  picture  qual¬ 
ity.  (See  figure  1 .) 

□  Improved  Basic 
Characteristics 

A  newly  developed  optical  structure  was 
introduced  in  these  devices.  This  struc¬ 
ture  adds  internal  lenses  above  the  photo 
shielding  film  to  the  Sony  Super  HAD 
CCD  technology.  (See  figure  2.)  By 
further  optimizing  the  lens  shape  using 
an  optical  simulator,  Sony  was  able  to 
increase  the  sensitivity  per  unit  area  by 
48%  over  Sony  1/4-inch  CCDs  with  the 
same  number  of  pixels,  thus  achieving 
a  level  of  300  mV, 

The  saturation  signal  level  was  also 
improved  by  65%  on  a  per  unit  area  ba¬ 
sis  by  optimizing  the  sensor  structure, 
achieving  values  of  600  mV  (NTSC) 
and  540  mV  (PAL).  Thus  these  devices 
achieve  excellent  performance  in  their 
imaging  c ha racteri sties  despite  being 
miniature  devices.  (See  table  2.) 


■  Further  miniaturization  and 
higher  picture  quality 

Due  to  the  development  of  the 
industry's  smallest  unit  pixel, 
these  devices  achieve  the  larg¬ 
est  number  of  pixels  in  a  1/6-inch 
optical  system  device. 

■  Improved  basic  characteristics 
The  sensitivity  per  unit  area  has 
been  increased  by  48%,  and  the 
saturation  signal  level  has  beer* 
increased  by  65%,  * 

■  Power  consumption  reduced  by 

38% 

*;  As  compered  (o  Sony  1/4-inch  devices  with  She 
same  number  of  pixels. 


O  Reduced  Power 
Consumption 

Current  consumption  was  reduced  by 
optimizing  the  output  circuits*  and  fur¬ 
thermore,  jmwer  consumption  was  re¬ 
duced  by  38%  over  Sony  1/4- inch 
CCDs  with  the  same  number  of  pixels 
by  reducing  the  capacitance  of  the  trans¬ 
fer  register.  These  devices  achieve  the 
lower  power  consumption  of  88  mW. 
including  the  driver  power  consump¬ 
tion,  (See  figure  3.) 

C  Device  Specifications 

The  specifications  of  these  devices. 
Including  the  number  of  pixels  and  the 
drive  specifications,  are  identical  to 
those  of  current  Sony  l /4-inch  CCDs 
with  the  same  number  of  pixels.  This 
means  that  current  system  [C  products 
can  be  used  without  modification.  (See 
table  L) 
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Micro  Lens  -  Sensor  Drawings 


Number  of  pixels  I  t  OK  pixels] 


Figure  1  Trends  in  CCD  Development 


1CX058CK  1CXQ65AK  ICX2Q8AK  1CX238AKE 

Produel  developed 
fnrlhB  relays 


Figure  3  Trends  in  Power  Consumption  in  38QK- 
pixel  CCDs 


Super  HAD  CCD  Sensor  Structure 


Figure  2  Sensor  Structure  Comparison 


Table  1  D e vice  S tru ctu re 


Item 

ICX230AKE 

ICX239AKE 

Optical  size 

116-inch  format 

1 /6-inch  format 

TV  formal 

NTSC 

PAL 

Transfer  method 

Interline  transfer 

InlerFine  transfer 

Total  number  of 
pixels 

811H*5G8V 

Approximately  41CK  pixels 

795H  x  596V 

Approximately  470K  pixels 

Number  of  effective 
pixels 

763H  x  494V 

Approximately  3S0K  pixels 

752H  x  502V 

Approximately  44 OK  pixels 

Chip  size 

3.30  mm  (H)  *  2.95  mm  (V) 

3.30  mm  {H)  x  2.95  mm  (V) 

Unit  cell  size 

3.200  tim  (H>  x  3.725  tim  (V) 

3.275  pm  (H)  x  3.150  pm  (V) 

Horizontal  drive 
frequency 

14.3132  MHz 

14.1875  MHz 

Package 

12  pin  SON  {Ceramic) 

9.25  mm  (H)  *  8.00  mm  (V)  x  2.30  mm  (t) 

12  pm  SON  (Ceramic) 

9.25  mm  (H)  *  8.00  mm  {V)  x  2.30  mm  (1) 

Table  2  Imaging  Characteristics 


Item 

ICX236AKE 

ICX239AKE 

Remarks 

Sensitivity  F5.6 

300  mV 

300  mV 

3200K.  706  mrf 

Saturation  signal 
level 

600  mV 

540  m  V 

Ta  =  60CP 

Smear  F5.6 

-86  dB 

-66  dR 

V/10  method 
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4.  Ray  Tracing  Code  and  Analysis 


A.  Ray  Tracing  Code  -  Introduction  to  Terminology 

IDA  developed  a  ray  tracing  code  (see  Appendix  A  of  this  document)  that  was  used 
to  analyze  optical  systems.  With  this  code,  we  calculated  the  MTF  and  PSF  to  get  a 
quantitative  comparison  of  the  resolution  of  different  optical  systems.  The  resolution  is 
determined  by  comparing  the  object  test  pattern  with  the  resulting  image  pattern  at  differ¬ 
ent  spatial  frequencies  (see  Figure  4-1). 

1.  MTF 

The  MTF  describes  the  loss  of  detail  in  the  imaging  process  as  the  convolution  of 
the  object  function  (i.e.,  the  ideal  image  function)  with  the  impulse  response  of  the 
imaging  system.  A  convolution  is  an  integral  that  expresses  the  amount  of  overlap  of  one 
function,  h,  as  it  is  shifted  over  another  function,/  It  “blends”  one  function  with  another. 
The  convolution  in  the  spatial  domain  is  a  multiplication  in  the  frequency  domain.  Thus, 
if/(x,y)  is  the  object  function  and  h(x, y)  is  the  impulse  response  of  the  imager,  then 

& g(x,y)]  =  <^[f(x,y)  x  h(x,y)],  (4-1) 

where  ^stands  for  Fourier  transform  and  leads  to 

Gfeq)  =  Ffeq)  *  H&q),  (4-2) 

where  x,  y  are  the  spatial  coordinates  and  /q  are  the  spatial  frequency  components  such 
that  if  X,  Y  are  the  spatial  periods  of  the  object  pattern,  then  £,  =1/X  and  q  =  1/Y. 
F  denotes  the  object  spectrum,  H(/q)  is  the  optical  transfer  function,  and  MTF  is  defined 
as  |H(£,,q)|. 

We  compared  the  MTF  of  optical  systems  obtained  with  our  code  and  the  published 
results  obtained  by  other  investigators.  In  addition  to  this  comparison,  we  also  compared 
our  code’s  results  with  analytical  solutions  for  simple  systems,  where  we  derive  an  ana¬ 
lytical  form  for  an  MTF  of  a  sine  function  with  a  Gaussian  PSF  and  calculate  the  MTF 
with  a  mathematical  program  following  the  analytical  closed-form  derivation.  This  com¬ 
parison  is  also  discussed  Appendix  A. 

2.  PSF 

The  PSF  describes  the  fuzzy  image  of  a  point  projected  by  the  optical  system  and  is 
represented  by  h(x, y)  for  real  systems.  For  ideal  systems  the  response  of  the  imager  is 
identical  to  the  object  function,  so  h(x, y)  =  </x,y)  or  the  delta  function.  The  PSF  is  related 
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to  the  MTF.  By  taking  the  Fourier  transform  in  two  dimensions  of  the  PSF  and  then 
taking  the  absolute  value  of  the  result,  we  obtain  the  MTF. 


Figure  4-1.  The  Spatial  Frequency  of  an  Object  Pattern 
With  the  Resulting  Image  Pattern  at  Different  Spatial  Frequencies 
Note  for  Figure  4-1:  This  figure  shows  spatial  frequency  variation  in  object  to  MTF  of  optical  system.  It 
shows  (a)  the  object  with  varying  special  frequencies,  (b)  the  variation  of  intensity  as  a  function  of  special 
frequency,  (c)  the  image  variation  of  intensity  as  a  result  of  an  imperfect  optical  system,  and  (d)  the  calcu¬ 
lated  MTF  of  the  system. 


B.  Briefing 

The  remainder  of  Chapter  4  presents  the  HARDI  Follow-Up  Meeting  (annotated) 
briefing  that  IDA  presented  to  DARPA/MTO  on  January  4,  201 1. 
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HARDI  Follow-Up  Meeting 
IDA  Support  to  HARDI 

Ray  Tracing  Analysis 

January  4,  2011 

Bohdan  Balko 
IDA 


For  this  meeting,  I  summarized  the  task  work  completed  before  October  5,  2010,  and  then 
described  a  ray  tracing  code  that  we  developed  and  verified  to  help  us  conduct  tradeoff 
studies. 
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Task  Objectives 

Evaluate  potential  benefits  of  introducing  curved  focal  plane  arrays 
(FPAs)  into  military  systems 

•  Simplified  and  substantially  improved  with  curved  FPAs 

•  Innovative  applications  (e.g.,  passive  three-dimensional 
(3-D)  imaging) 

•  Gradient  Index  of  Refraction  (GRIN)  lenses  with  curved  FPAs 


The  question  is,  how  can  we  improve  military  systems  by  reducing  the  size,  weight, 
complexity,  and,  therefore,  the  cost  of  existing  systems?  What  are  some  innovative  ways  to 
solve  problems  previously  not  possible  with  plane  focal  plane  arrays  (FPAs)?  How  can  we 
combine  Gradient  Index  of  Refraction  (GRIN)  lenses  with  curved  FPAs  in  an  optimized 
way? 
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Suggested  Military  Applications 

•  Small  robots  used  in  the  Advanced  Mine  Detection  System 
(AMDS)  program 

•  Robots  used  in  the  Explosive  Ordnance  Disposal  (EOD) 
program 

•  Miniature  unmanned  aerial  vehicle  (UAV) 

•  Small  surveillance  camera 

•  Spectral  bandwidth  opening  up  new  applications 


We  found  military  applications  and  presented  them  at  the  first  Hemispheric  Array  Detector 
for  Imaging  (HARDI)  follow-up  meeting  that  I  attended  in  July  2010.  The  sponsor  selected 
small  robots  used  in  mine  detection  and  instructed  IDA  to  leam  about  the  characteristics  of 
the  cameras  used  on  the  robots.  Specifically,  the  HARDI  team  wanted  to  know  the  pixel  size 
or  pitch,  the  number  of  pixels  in  the  FPA,  and  the  field  of  view  (FOV),  along  with  other 
parameters  (e.g.,  dimensions). 
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Robots  and  Their  Cameras 

•  i Robot  [PackBot,  Warrior) 

Contact:  Orin  Hoffman 

•  QinetiQ  North  America  (formerly  Foster  Miller)  (Talon) 
Contact:  Peter  Wells,  Senior  Engineer 


The  robot  manufacturers  and  their  representatives  listed  in  this  slide  provided  the 
information. 
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Cameras  Used  on  Robots 


PC168-IR  on  Talon 

FOV:  80  deg 
Pixel  size:  2.85  pm 


Zoom  Camera  on  Talon 

Pixels:  680,000 
FOV:  42,2  deg 


i  Robot  ( PackBot ,  M/arr/or) 

Sony  FCB-EX780R 
Pixel  size:  3.2  pm  x  3.725  pm 
Pixels:  680,000 

FOV:  45  deg  (wide  end),  2,0  deg  (tele  end) 
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The  three  cameras  about  which  I  was  able  to  obtain  information  are  shown  in  this  slide. 
More  detailed  specification  sheets  were  appended  to  the  original  briefing  presented  to  the 
HARDI  group  in  October  2010  (see  Chapter  3  of  this  document). 
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Other  Activity 

•  I  attended  QinetiQ  robot  demonstration  at  Quantico  Marine 
Base  on  August  3,  2010.  My  interest  in  attending  the  demo 
was  to  see  the  robots  in  action  and  obtain  information  about 
the  cameras  used  on  the  robots 

•  One  robot  was  equipped  with  two  cameras  to  produce  3-D 
imaging.  I  tried  to  control  the  robot's  grips  with  and  without 
3-D.  It  was  clearly  beneficial  to  have  3-D.  I  mention  this 
because  there  may  be  an  innovative  way  to  provide  passive, 
3-D  target  imaging  with  curved  F PA  and  one  camera 
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Future  Proposed  Work  (10/5/2010) 

•  Conduct  a  trade  study  of  pixel  size  and  count  vs.  camera 
characteristics  to  obtain  optimum  conditions  for  various 
robot  camera  applications 

-  To  prepare  for  this  trade  study,  we  developed  a  ray  tracing 
code  and  investigated  a  well-known  recent  publication  on 
curved  detectors 

•  Continue  to  obtain  information  about  cameras  used  on 
robots  in  AMDS  and  EOD  programs  to  get  the  latest  specs 

•  Compare  pixel  material  developed  by  the  Hemispheric  Array 
Detector  for  Imaging  (HARDI)  program  with  material  in 
current  cameras.  Look  at  the  result  if  HARDI  material  was 
used  in  a  flat  FPA  for  a  direct  comparison 


At  the  October  5th  HARDI  follow-up  meeting,  I  proposed  to  conduct  a  tradeoff  study  to 
determine  the  optimum  parameters  for  the  HARDI  systems  to  be  designed  for  various 
applications. 

To  do  this,  we  needed  a  tool — computer  code — to  study  the  effects  of  various  parameters  that 
affect  the  system  design. 

We  studied  the  characteristics  of  spherical  FPA  published  by  the  Stanford  Group  and  intended 
to  compare  the  results  of  the  IDA  code  using  the  same  parameters  and  geometry. 
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Analysis  of  Stanford  Paper  on  Curved  FPAs 

Rim,  S-B.,  B.  Catrysse,  R.  Dinyari,  K.  Huang,  and  P.  Peumans.  2008. 
"The  Optical  Advantages  of  Curved  Focal  Plane  Arrays."  Optics 
Express  16  (7):  4965-4971. 
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Stanford  Publication 

•  The  authors  investigate  the  potential  advantages  to  optical  imagers  of  replacing  the  standard  planar 
imaging  surface  with  a  hemispherical  detector 

*  The  authors  assert  that  with  spherical  symmetry  [with  concentric  optics  and  detector),  any  object 
point  can  be  considered  "on-axis" 

-  This  would  remove  the  distinction  between  tangential  and  sagittal  rays 

-  it  would  also  eliminate  four  out  of  five  Seidel  aberrations  [coma,  astigmatism,  field  curvature,  and  distortion) 

*  System  is  ostensibly  designed  for  wide  field-of-view  [FOV)  applications 

-  Sym m et ry  becomes  more  adva n tageo us  as  F 0 V  i n c reases 

-  Wit  h  th  e  sma  1 1  f  oca  I  len  gt  hs  be  i  ng  co  n  sidered  he  re  [~6  mm),  a  ngu  I  ar  resol  utio  n  wi II  n e ve r  a p p  roach  t  he 
capabilities  of  a  modern  narrow  FOV  imager 

•  Authors  have  modeled  the  proposed  system  and  compared  it  to  two  other  systems.  At  f/3,5,  the 
proposed  system  seems  to  outperform  the  others: 


Plano-convex:  Cooke  Triplet:  Concentric  Spheres: 


The  figures  at  the  bottom  of  this  slide  give  the  MTFs  calculated  by  the  Stanford  group  for  three 
different  lens  systems.  The  view  on  the  right  shows  the  MTF  for  the  ball  lens.  An  interesting 
observation  is  a  comparison  between  the  diagram  for  the  ball  lens  and  the  two  more  traditional 
systems.  The  left  and  middle  figures  appear  to  present  rays  that  can  pass  anywhere  through  the 
lens  systems.  However,  in  the  figure  on  the  right,  all  the  displayed  rays  traverse  the  lens  fairly 
close  to  its  center,  as  if  an  aperture  were  located  in  the  center  of  the  ball.  The  paper  does  not 
mention  an  aperture  in  the  center  of  the  nor  does  it  discuss  the  different  constraints  imposed  on 
the  systems. 
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Code  Development  at  IDA 

•  To  study  the  effects  of  the  various  system  parameters  and 
investigate  the  tradeoffs,  we  developed  a  ray  tracing  code  for 
calculating  the  modulation  transfer  function  (MTF)  ofsystems- 
of-interest,  along  with  their  characteristics  (MATLAB) 

•  We  checked  the  code  by  comparing  it  with  published  results 
{i.e.,  Rim,  S-B.,  et  al.,  "The  Optical  Advantages  of  Curved  Focal 
Plane  Arrays") 

•  Two  system  geometries  used  with  different  aperture  positions 

•  We  also  obtained  an  analytical  (closed-form)  MTF  result  for  a 
simple  system  and  studied  the  result  for  various  spatial 
frequencies  and  point  spread  functions  (PSFs)  using 
Mathematica  and  Mathcad 

•  Satisfied  that  the  MATLAB  code  gives  reasonable  results,  we 
are  ready  to  conduct  tradeoff  studies 

12 


This  slide  describes  the  IDA  code  development  and  verification  procedures. 
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IDA  Ray  Tracing  Model  Characteristics 

Characteristics  of  base  model: 

•  Diffraction  is  ignored 

•  The  lens  is  a  sphere,  characterized  by  a  single  real  index  of  refraction 

•  Computations  are  performed  in  one  plane  (meridional) 

•  Meridional  PSFs  and  MTFs  are  calculated 

We  can  include  diffraction,  dispersion,  skew  rays,  multi-dimensional 
PSFs  and  MTFs,  and  more  general  optical  materials  {GRIN  lenses)  and 
configurations  in  subsequent  versions.  We  can  also  modify  the  FPA 
surface  to  study  more  general  shapes  {i.e.,  parabolic  FPAs) 


This  slide  describes  technical  characteristics  and  assumptions  of  the  IDA  model. 
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Simple  Example  of  Results  Form  IDA  Model 


Aperture  width  -  180  deg 
Object  Angle  =  on  axis 
Object  distance  »  Lens  radius 
Index  of  refraction:  nlENS  =  1.5 
Lens  radius:  rLENS  =  4,0  mm 
Detector  radius:  rDETecr0K  =  6,0  mm 
(i,e,,  detector  is  at  focal  position) 


Diagram  of  Optical  System  and  Detector 


Poinl  Spread  Function .  Optical  Aberrations 


Modulation  Transler  Function,  Optical  Aberrations 


This  example  shows  some  of  the  output  that  we  can  produce  with  our  code.  The  upper  left 
box  shows  the  parameters  used  in  the  calculation.  First,  we  produce  the  ray  traces  shown  in 
lower  left  diagram.  Although  we  use  many  more  rays  than  shown,  we  only  display  these  few 
selected  ones  to  prevent  confusion.  The  upper  right  diagram  is  the  PSF,  and  the  lower  right 
diagram  is  the  aberration  MTF  (blue  line)  and  the  detector  pixel  MTF  (blue  line). 
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This  slide  shows  two  different  geometries  used  in  our  calculations.  The  difference  between 
these  geometries  is  the  position  of  the  aperture.  Geometry  I  uses  an  aperture  at  the  ray  entry 
to  the  lens.  Geometry  II  places  the  aperture  in  the  plane  crossing  the  center  of  the  ball  lens. 
In  Geometry  II,  all  the  rays  are  effectively  close  to  axial.  For  Geometry  I  calculations,  we 
used  an  f#  given  by  focal  length  divided  by  aperture  diameter.  For  Geometry  II,  we  have 
redefined  f#  based  on  the  slope  of  the  exit  cone  only. 
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Defining  F#  for  Thick  and  Thin  Lenses 


Appropriate  Definition 

for  Ball  Lens 
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This  slide  shows  the  geometric  parameters  traditionally  used  to  calculate  F#  for  a  thin  lens 
system  (above)  and  what  we  feel  is  appropriate  for  the  ball  lens  (below). 

This  is  a  crucial  criterion  to  match  appropriately  because  the  spatial  extent  of  aberration 
PSF  is  proportional  to  1  (in  the  small-angle  limit). 

(F#)3 
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IDA  Calculation:  Geometry  I 


*  For  the  conditions  specified  in  the  Stanford  paper,  the  IDA  model  for  Geometry  I  gives 
reasonable  results,  although  somewhat  less  sanguine  than  those  presented  earlier: 


R#  -  3,5;  identical  to  Stanford 

Lens  radius  -  4  mm:  Identical  to  Stanford 

nLEws "  1-S168  (nd  for  6K-7  glass);  Stanford  used  a 

weighted  sum  of  the  response  at  Fraunhofer  lines  C 

(H),  F  (H),  and  d  (He) 

nc  =  1.51432,  n f  =  1.5223 8,  nd  =  1,51680 

Detector  radius  (and  position)  =  5.87  mm:  Identical  to 

Stanford,  Note  that  this  is  the  focal  length  of  the  lens 

(at  nd) 


Diagram  of  Optical  System  and  Detector 


Mud  illation  Transfer  Fun^tioni,  Optical  Aberrations 


Spatial  Frequency  ^cycJmm] 


The  Stanford  paper  shows  an  MTF  of  about  0.6  (-0.75  of  the  diffraction  limit)  at 
68  cyc/mm.  Our  computation  yielded  an  MTF  of  0.42  at  68  cyc/mm. 
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Geometry  I:  Off-Axis  Behavior 

•  We  have  only  considered  rays  at  the  center  of  the  FOV.  What  happens  as  we  move 
toward  the  edges  of  the  array? 

-  The  rays  strike  the  aperture  at  a  glancing  angle  and  are  asymmetrically  deviated 

-  Coma  and  astigmatism  will  be  large 

-  In  this  case,  aberration  would  be  reduced  by  moving  the  detector  much  closer  to  the  lens 


Modulation  Transfer  Function,  Oplical  Aberrations 


The  symmetry  has  been  broken  by  the  aperture.  Will  this  optical  system  really  work? 
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Geometry  I  gives  an  awful  MTF  result  off  axis.  The  straight  rays  are  blocked,  and  the  ones 
that  are  permitted  are  sharply  bent,  leading  to  significant  aberration. 
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IDA  Calculation;  Geometry  II 

*  For  the  conditions  specified  in  the  Stanford  paper,  the  IDA  model  gives  comparable  but 
somewhat  less  optimistic  results: 


F#  =  3*5:  Identical  to  Stanford 

Lens  radius  -  4  mm:  identical  to  Stanford 

nlEM5=  1.5168  |nd  for  BK-7  glass):  Stanford  used  a 

weighted  sum  of  the  response  at  Fraunhofer  lines  C 

(H),  F  (H),  and  d  (He) 

nc  -  1.51432,  nF  -  1,52238,  nd  - 151680 

Detector  radius  (and  position)  =  5.87  mm:  Identical  to 

Stanford.  Note  that  this  is  the  focal  length  of  the  lens 

(at  nd) 


Diagram  of  Optical  System  and  Del  actor 


Point  Spread  F  Line  I  ion.  Qplical  Abcrrallons 


Modulation*  Transfer  Funeikm,  Optical  Aberrations 


The  Stanford  paper  shows  an  MTF  of  about  0.6  (-0.75  of  the  diffraction  limit)  at 
68  cyc/mm.  Our  computation  yielded  an  MTF  of  0.47  at  68  cyc/mm. 

Geometry  II  does  slightly  better  than  Geometry  I,  although  the  MTF  is  still  not  as  high  as 
that  shown  in  the  Stanford  paper. 

The  difference  between  Geometry  I  and  Geometry  II  on  axis  is  only  because  of  the 
differences  in  f#  definition. 
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Geometry  II:  Off-Axis  Behavior 

*  What  happens  as  we  move  toward  the  edges  of  the  array? 

-  The  projection  of  the  hole  appears  smaller  at  45  deg 

-  PSF  narrows,  so  MTF  improves  (at  the  cost  of  a  decrease  in  signal  strength) 


Diagram  of  Optical  System  and  Detector 


Modulation  Transfer  Function,  Optical  Aberrations 


This  is  an  enormous  improvement  over  Geometry  I  ! 

In  fact.  It  is  even  better  than  the  on-axis  performance  for  either  geometry 

20 


Geometry  II  preserves  (and  even  improves)  the  MTF  off  axis. 
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Adjustment  of  Detector  Position 

Geometry  II  is  clearly  the  better  way  forward.  However,  two  questions  remain: 

*  Why  are  the  published  MTFs,  which  presumably  include  diffraction,  higher  than  ours? 

*  Are  there  any  simple  adjustments  we  can  make  to  increase  our  MTFs? 

Also,  the  focal  position  is  a  sub-optimal  location  for  the  detector 

*  In  a  spherically  aberrated  converging  optical  system,  some  benefit  is  gained  in  moving  the 
detector  slightly  closer  to  the  lens  than  the  focal  position.  The  optimal  point  trades  focal 
precision  of  rays  nearest  the  center  of  the  entrance  pupil  for  minimization  of  total  blur 
spot  area  (circle  of  least  confusion) 


The  circle  of  least 
confusion  occurs 
where  the  marginal 
rays  exit  the  caustic 
envelope.  This  occurs 
at  the  location  XLC  in 
the  figure* 


*  Figure  from  Eugene  Hecht,  Optics ,  3rd  ed.  (Reading,  MA:  Addison  Wesley  Publishing  Company  1997)  21 


The  circle  of  least  confusion  is  the  position  where  the  image  plane  (detector)  should  be  placed 
to  minimize  the  overall  size  of  the  aberration  blur  spot.  For  a  converging  lens  with  positive 
longitudinal  spherical  aberration  (LSA),  as  is  the  case  here,  this  position  will  be  located  in 
front  of  the  focal  position.  Assuming  axial  symmetry,  the  circle  of  least  confusion  is  defined  as 
the  locus  of  points  where  a  marginal  ray  exits  the  caustic  ray  bundle  (envelope). 


4-21 


Optimal  Detector  Position 


MTF  improves  significantly 
when  the  detector  is  placed  at 
the  position  of  least  confusion 

This  placement  improves  the 
overall  MTF  and  greatly  reduces 
its  dependence  on  object 
position 


Optical  AJmhi attar,  MTF,  DatocHsr  tC  FkjI  PmMton  (5.S7  nvn| 


Ablation  HTF*.  DHST^nt  Dd'bcEqt  Ft** 


Qplic.il  uiiarr.ilioii  MT  F,  OolfilH  -g  Poitlun  nl  L*ai1  CDfifbshHi  |  S  U  j  mm| 


By  changing  the  position  of  the  FPA,  we  can  get  a  variety  of  MTF  results.  For  the  focal 
position  at  5.87  mm,  which  was  the  position  used  by  the  Stanford  Group,  we  compute  about 
0.47  for  the  MTF  at  68  cyc/mm.  The  position  at  5.83  mm  (circle  of  least  confusion)  gives 
0.8  for  the  MTF  at  68  cyc/mm — a  much  better  result. 

The  two  charts  at  the  bottom  of  this  slide  show  how  the  dependence  of  MTF  upon  object 
angle  is  also  diminished  by  moving  to  the  circle  of  least  confusion. 
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Possible  Enhancements 

The  following  considerations,  while  not  mentioned  by  the  authors, 

might  also  be  useful: 

1.  Could  the  lens  be  coated  with  a  material  that  would  transmit  rays 
near  normal  incidence  but  would  strongly  reflect  radiation  at 
other  angles? 

2.  Could  a  spherically  symmetric  GRIN  lens  allow  a  wider  aperture? 

3.  Could  the  detector  position  (and  size)  be  optimized  further,  in 
combination  with  Options  1  and  2? 
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Conclusions 

We  have  demonstrated  that  the  proposed  optical  system  could 
lead  to  consistent  resolution  across  a  wide  FOV 

The  ideal  location  for  the  detector  is  not  the  focal  position  but 
the  position  of  least  confusion 

The  requirement  for  an  aperture  (for  a  reasonable  blur  spot) 
breaks  the  spherical  symmetry  alluded  to  in  the  Stanford  paper, 
but  the  effect  on  MTF  is  small  at  the  position  of  least  confusion 
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Symmetry  Restored:  Geometry  1 

Symmetry  could  be  restored  by  opening  up  the  aperture  as  widely  as  possible ... 

At  f/0.75,  all  straight  rays  are  permitted  to  pass  through  the  system.  Only  at  the  very 
edges  will  the  blur  spot  and  MTF  change  appreciably.  However,  the  resulting 
aberrations  are  awful 

„  On-axis  MTF 

9^*  *  "  [ 

'L _ ^  v _ _J  I- 

ir4=ar**gll 

1 

.  tfl  Opfccdl  aTl1iMTi  Mrf 

45  deg  off-axis  MTF 

.A.  1 

1 

. 

— - - 

r  *  Se 

We  initially  investigated  the  authors’  symmetry  claims  with  Geometry  I  and  found  results 
that  disagreed  strongly.  The  reason  for  the  discrepancy  was  that  the  insertion  of  an  aperture 
had  broken  the  spherical  symmetry  that  was  purported  to  ensure  that  the  image  resolution 
would  not  deteriorate  as  the  object  moved  off  axis  toward  the  edges  of  the  FOV. 

Here,  we  note  that  by  eliminating  the  Geometry  I  aperture,  we  can  preserve  the  spherical 
symmetry.  However,  the  aberrations  now  dominate,  and  the  results  are  poor. 
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By  changing  the  position  of  the  FPA,  we  can  get  a  variety  of  MTF  results.  The  position  at  5.87 
mm,  which  was  used  by  the  Stanford  Group,  gives  about  0.47  for  the  MTF  at 
68  cyc/mm.  The  position  at  5.83  mm  (circle  of  least  confusion)  gives  0.8  for  the  MTF  at  68 
cyc/mm — a  much  better  result. 
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Circle  of  Least  Confusion 

Light  rays  lying  in  the  tangential  and  sagittal  planes  are  refracted  differently  and  both 
sets  of  rays  intersect  the  chief  ray  at  different  image  points,  termed  the  tangential 
line  image  (tangential  focal  plane)  and  the  sagittal  line  image  (sagittal  focal  plane). 
These  rays  fail  to  produce  a  focused  image  point;  rather,  they  produce  a  series  of 
elongated  images  ranging  from  linear  to  elliptical,  depending  upon  the  position 
within  the  optical  train.  In  a  zone  known  as  the  circle  of  least  confusion  (positioned 
between  the  tangential  line  image  and  the  sagittal  line  image),  the  major  and  minor 
axes  of  the  ellipse  are  equal  and  the  image  approaches  a  circular  geometry 


Source:  h ttp :// m f c ro. m ag net .tsu.edu/prinne r/j a va/a be rrat i on s/ astigmat ism/ index. h tmi 


The  circle  of  least  confusion  is  the  smallest  image  (blur  patch  of  light)  of  point  source  that 
optical  system  can  project.  As  lens  moves  in  and  out  of  focus,  patch  of  light  grows  and 
decreases  in  size.  The  smallest  size  is  circle  of  least  confusion,  and  the  optical  system  is  said  to 
be  in  focus.  Note:  due  to  the  wave  nature  of  light  and  depending  on  lens  design,  the  circle  of 
least  confusion  does  not  necessarily  give  best  definition  (Source: 
http://www.idigitalphoto.com/dictionary/circle  of  least  confusion.) 
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Appendix  A. 
Ray  Tracing  Codes 


Analytical  Modulation  Transfer  Function  (MTF)  Calculation 

In  the  process  of  developing  the  main  ray  tracing  code  for  our  general  use,  we 
decided  that  we  needed  a  simple  closed  form  result  to  check  our  numerical  results.  This 
would  provide  assurance  that  the  code  was  dependable — at  least  for  simple  functions. 

This  is  the  closed-form  solution  used  to  calculate  the  MTF  for  simple  functions  (i.e., 
sine  with  (1)  varying  spatial  frequency  as  the  object  and  (2)  Gaussian  of  various  widths 
as  the  point  spread  function  (PSF)).  This  result  was  used  to  check  the  MATLAB  code 
RAYMTF  and  gain  confidence  in  the  RAYMTF  code’s  MTF  results  for  more  complex 
systems. 

We  start  with  two  functions:  a  sine  wave  with  frequency  f,  Sin(fx),  and  a  Gaussian 
with  parameter,  e~ax 2  ,  as  the  impulse  response.  The  parameter  “a”  is  related  to  the  typi¬ 
cal  width  of  the  Gaussian  by  the  relationship  =  —  . 

The  first  step  is  to  do  the  convolution  of  the  two  functions: 

j  dx  Sin(fx)e~a^x~t^2 .  (A-l) 


To  facilitate  the  calculation,  the  exponential  description  of  the  sine  is  used: 

j+ifx-a(x-t)2  _  0-ifx-a(x-t )2 

J 

_  —  J  dx  ^g-ax2-(2at-if)x-at2  _  g-ax2-(2at+if)x-at (A-2) 


2-J  dx  (e+ifx  -  e-Vx)e-a(.x-t)2  =  2_  (e+i/*-a(*-Oz  -  e-ifx-a(x-t)^ 


Focusing  on  the  first  term,  we  can  solve  it  by  completing  the  square  to  get  a  shifted  Gaus¬ 
sian  and  a  constant  exponential: 


-ax2  -  (2 at  —  i f)x  —  at2  =  -a  [x2  —  2  (t  +  i x  +  t2 j 

=  —  a  |x2  —  2  (t  +  i  +  (t  +  i-^j  -  (t  +  /A)  +  t2 

=  -a[(x-(t  +  i£))  ~{t2  +  l-^-{^)  +  t2] 

=  -a(*-[t  +  i£])  +(+i/t-A). 


(A-3) 


Going  back  to  the  integral,  we  get 


A-l 


(A-5) 


r  dx  a(x-[t+i£])2  +ift-f2 /4a 
J  2  i 

The  only  difference  between  the  first  and  second  terms  is  the  sign  of  the  complex 
term,  ift,  so  by  making  the  substitution  in  the  previous  calculation  of 

+ift  ->  —ift ,  (A-6) 


we  get  for  the  second  term 

f  ^ e-a(x-[t-fA])2  e_i/t_/2/4a  (A_7) 

The  two  terms  can  now  be  put  together  (the  Gaussian  integral  calculated  and  con¬ 
stant  terms  rearranged)  to  get 


1 — 
n 
-( 

\ a 


-2  i 


(A-8) 


To  get  the  MTF,  we  divide  by  the  original  function,  Sin  (ft),  to  get 


MTF(f,a)  =  J^e~/2/4a- 


(A-9) 


a  = 


To  relate  back  to  the  original  width  of  the  Gaussian,  we  use  the  relationship 

4<J2 ' 


MTF(f,a)  =  J^e~(/cr)2 


(A- 10) 


(This  space  intentionally  left  blank.) 
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MTFCHECK1 


This  is  the  “Mathematica”  code  used  to  graph  the  analytical  MTF  results. 

F[x_,  f_]  :=  Sin[fx  /  (2  Pi) ]  ; 

G[x_,  a_]  :=  Sqrt [Pi  /  a]  Exp  [-a  xA 2] ; 

F [x,  f] 


Sin 


G[x,  a] 


Convolve  [Sin  [f  x  /  (2  Pi) ]  ,  Sqrt  [Pi  /  a]  Exp[-axA2]  ,  x,  y] 

$Aborted 

F[x,  f] 


f  =  2;  a  =  1; 

F[x,  f] 


Convolve  [  x,  1/  (xA2  +  1)  ,  x,  y] 

Try 

Convolve  [Sin  [  ff  x]  ,  E  A  (-axA2)  ,  x,  y] 


1 

- ie'< 

2 


--ff  (ff+4  i  y) 


(- 1  +  e2 1  ff y)  VtT 


--ff  (ff+4  ±  y) 


(-1  +  e2iff y) 


ComplexExpand [h [y ] ] 

1  _ _  1  _ 

—  e  *  Vtt  Sin  [f  f  y] - e~  *  Vrr  Cos  [2  f  f  y]  Sin  [ff  y] 


1 


1 

+  —  e 


2 


2 


2 


1 

i  —  e* 

2 


2 


2 
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ComplexExpand [h [y ]  /  Sin [f f  y]  ] 


ff2  _  ff2  _ 

e"T  V7T  Sin[ff  y]2  e"“  VtT  Cos  [2  ff  y]  Sin[ff  y]2 


- 1  +  Cos  [2  ff  y] 


-1  +  Cos  [2  ff  y] 


/ — 

e  «  Vtt  Cos  [ff  y]  Sin  [ff  y]  Sin  [2  ff  y] 
- 1  +  Cos  [2  ff  y] 


e  *  VtT  Cos[ffy]  Sin[ffy] 
-1  +  Cos  [2  ff  y] 


ft2  _  ff2  _ 

e“  «  Vtt  Cos  [f  f  y]  Cos  [2  ff  y]  Sin  [ff  y]  e~~  Vtt  Sin  [f f  y] 2  Sin  [2  f f  y] 


-1  +  Cos  [2  ff  y] 


-1  +  Cos  [2  ff  y] 


Plot [  (h [y] )  /.  (ff  ->  10)  ,  {y,  0 ,  10}] 

Plot [  (Sin [ffy])  /.  (ff->10),  (y ,  0,  10}] 


ExpToTrig  [h  [y]  /  Sin  [  f f  y]  ] 

- i  VtT  Csc  [f f  y]  (- 1  +  Cos  [2  f f  y]  +  i  Sin  [2  f f  y] )  (cosh  [  —  f f  (f f  +  4  i  y)  1  -  Sinh  [  —  ff  (f f  +  4  i  y) 

2  {  L  4  J  L  4 


a  =  1 

1 

Convolve  [Sin  [f  x]  ,  EA(-bxA2),  x,  y.  Assumptions  ->  b  >  0  ]  /  (Sin  [f  y] ) 

i  e~  '  <»  VtT  Csc  [f  y]  2  +  2  e2 1  f  ^  t  Erf  [  y|>i  Erf  i  [  ]  | 

4  Vb~ 
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ExpToTrig [ - 


(f»4  Iby)  _  |  . 

ie  4b  V7r”csc[fy]  |-2  +  2  e2if  Y  +  Erf  +  i  Erfi  f*2^y_  j  j 

4  Vb” 


4  Vb” 

i  VtT  Csc  [f  y] 


-2  +  2  Cos  [2  f  y]  +  Erf  [ - - - +  Vb~  yl  +  i  Erfi  [ - +  i  VtT  yl  +  2  i  Sin [2  f  y] 

2  Vb~  2  Vb” 


r  f  (f  +  4  1  b  y)  1 

r  f  (f  +  4iby) 

Cosh  - 

-  Sinh 

L  4.b  J 

L  4  b  J  j 

Simplify [%,  f  >  0  &&  b  >  0  &&  Element  [b,  Reals]] 

1 


4  Vb~ 


i  VtT  Csc  [f  y] 


r-if+2by,  rf+2iby1 

2  +  2  Cos  [2  f  y]  +  Erf  [ - J  +  i  Erf i  [ - J  +  2  i  Sin  [2  f  y] 


2  Vb~ 


2  Vb” 


j  rf  (f+4iby)1 

rf  (f +  4iby) 

Cosh 

Sinh 

l  L  4b  J 

L  4  b  J  j 

Simplify  [ 
1 


4  Vb” 


£i  n 

-e"2o  j  — 


ie  4b  Vtt"  Csc[fy] 


r-if +  2by,  rf  +  2iby, 

-  2  +  2  e2  Af  y  +  Erf  [ - J  +  i  Erf i \ - 1 


2  Vb” 


2  Vb” 


/-  b-*  5 


Plot 


rf  1  f(f*4*by)  y—  f  r-if +  2by,  rf +  2iby,] 

- ie  4b  Vtt  Csc  [f  y]  -2  +  2  e2  Af  y  +  Erf  [ - 1  +  i  Erf  i  | - 1 

LV4 vs-  l  L  2vr  J  L  2vr  JJ 

{f,  0,  100} ,  PlotRange Allj 


/.  b  ->  1, 


1.0 

0.5 

i_i_ i_ i_ i_ i_ i_ i , i 

•  >  •  .  1  .  .  .  .  1 

•o 

© 

1 

o 

7 

0.5  1.0 

-0.5 

-1.0 

MFT1  [f  ,  b  ]  :  =  e~4b 
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Plot  [Evaluate  [Table  [MFTl[f,  b]  ,  (b,  1,  100,  20}]],  (f,  0,  20),  PlotRange -*  All] 


Plot[Sqrt[Pi  /  a]  EA(-axA2),  (x,  -5,  5},  PlotRange -»  All] 


Convolve [HeavisidePi [x] ,  Sin[x],  x,  y] 

Convolve [HeavisidePi [x] ,  Sin[x],  x,  y] 


(This  space  intentionally  left  blank.) 
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MTFCHECK2 

This  is  the  “MATHCAD”  code  used  to  graph  the  analytical  MTF  results. 


Af  :=  0.1 
T  :=  fO  +  Af-j 


fO  :=  0.1 


j  :=  1,2..  700 
a  :=  0.007 


MTF. 

J 


a 


Act  :=  0.01 


ctO  :=  0.0001 


(1^:=  ctO  +  Acj-k 


MTF2.  , 

J»k 


k:=  0,1..  10 

sin(fixi) 


Ax:=  0.01 

1  :=  0,1..  100 

1*2..  10 

X|  :=  Ax-1 
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cO  :=  O.OOOOl  Aa  :=  0.1  aO  :=  0.0001 
Ac  :=  0.002 

c  :=  cO  +  Ac-m  a  :=  aO  +  Aam 

m  m 

Xj  =  0.01 


G  :=(l+c  .^)e'(0'5‘X|^2'(Cm)2] 
l,m  \  m  v  / 


a5  =  0.5 


CT5  =  0.05 


(This  space  intentionally  left  blank.) 


A-8 


RAYMTF 


This  is  the  MATLAB  code  used  to  calculate  the  ray  tracing  results  and  the  MTFs 
and  PSFs  shown  in  Chapter  4. 

*  rek  20110120 
%to  do:  flat  fpa 
%to  do:  thin  lens 
%to  do:  compound  optics 

%to  do:  allow  multiple  wavelengths...  (dispersion) 

%  1)  INPUT  PARAMETERS: 

%  la)  SCENARIO: 

numObjectPoints  =  1; 

objectPoints_def inedbyPolarCoords  =  true; 

objectAngles  =  180*pi/180;  %  radians  %  (if  using  POLAR  Coords) 

%  measured  from  +x  axis  (row  vector  of  Obj Angles,  one  per*^ 

Ob j Point) 

objectDistance  =  1000;  %  meters  %  (if  using  POLAR  Coords) 

%  if  very  big  compared  to  sensor  size,  incoming^ 

rays  will  be  parallel 

objectPoints _ CartesianDefault  =  [-le3;0];  %  meters  %  (if  not  using*' 

POLAR  Coords) 

%  lb)  SENSOR  PARAMETERS: 

%  lbl)  optical  elements: 
lensRadius  =  0.004;  %  meters 

lensCenter  =  [ 0 ; 0 ] ;  %  meters  %  Cartesian  coordinates  of  lens  center 

nLens  =  1.5168;  %index  of  refraction  n_D3  (i.e.  at  yellow  Helium  line  at  587.6  nm) 

%  note  that  Stanford  weighted  their  design,  analysis  and  results*' 

for  n_C,  n_D3,  and  n_F: 

%  weights  =  {1,2,1} 

%  according  to  http://refractiveindex.info/?*' 

group=GLASSES&material=BK7 

%  n_C  =  1.51432  @  656.3  nm 
%  n_D3  (=  n_d)  =  1.5168  @  587.6  nm 
%  n_F  =  1.52238  @  486.1  nm 

nAir  =  1;  %index  of  refraction  for  surrounding  medium 
%  lb2)  aperture  stops: 


useFnumberinput_toCalculateFrontsideAperture  =  true; 

front  of  lens 

%  Geometry 

I: 

aperture  at ^ 

frontsideAp _ f Number  =0; 

f Number  input) 

%  Geometry 

I 

(if  using*' 

frontsideApertureRadius _ default  =  lensRadius; 

using  f Number  input) 

%  Geometry 

I 

(if  not*' 

useFnumberinput_toCalculateInsideAperture  =  true; 

inside  lens 

%  Geometry 

II: 

:  aperture*' 

insideAp _ fNumber  =  3.5; 

f Number  input) 

%  Geometry 

II 

(if  using*' 

insideApertureRadius _ default  =  lensRadius;  %  Geometry  II  (if  not*' 

using  f Number  input) 

%  lb3)  detector: 

fpaCenter  =  [ 0 ; 0 ] ;  %  Cartesian  coordinates  of  fpa  center 

putFPAatFocalPosition  =  false; 

designobjectDistance  =  Inf;  %ObjDist  which  determines  FPA  position  (use  lens*' 
maker's  eqn) . 

putFPAatCircleOfLeastConfusion  =  false; 
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hardwireFPAposition  =  not (putFPAatFocalPosition  II  putFPAatCircleOfLeastConfusion) ; 

detectorRadius  =  0.00583;  %  meters  %  (only  used  if  hardwiring  FPA/ 

position) 

%  lc)  PRECISION  of  Calculation  (and  Charts) : 

numRays_perObjectPoint  =  10000;  '^number  of  evenly  spaced  meridional  rays  to  draw  from*' 
the  object  to  the  lens 

nPSFsampleBins  =  65536;  *the  number  of  bins  with  which  to  sample  the  optical  point / 
spread  function 

maxPlot  =  105;  %limit  for  number  of  rays  to  plot 

pixWidth  =  pi/1000;  %angular  size  of  pixel  used  for  computing  detector  MTF  (radians) 

%  End  input  parameters 


%  2)  GET  CONFIG  INFO  (Get  All  Configuration  Info  from  Input  Parameters): 

%  2a)  object  Points: 
if  ( ob j ectPoints_def inedbyPolarCoords ) 

objectPoints  =  objectDistance* (cos (object Angles) ; sin (objectAngles) ] ;  coordinates  of*' 

point  sources 
else 

objectPoints  «  objectPoints _ CartesianDefault ; 

end 

%  2b)  focal  position: 

f ocalRadius_ordl  =  1/ ( (2* (nLens-1) ) / (nLens*lensRadius) -1/designobjectDistance) ; 

%lens  makers  equation  (focuses  objects  at  "designobjectDistance" ) / 

FIRST  ORDER  ! 

%  2cl)  front  aperture: 

if  (useFnumberinput_toCalculateFrontsideAperture) 

%  see  comments  on  ’ useFnumberinput_toCalculateInsideAperture •  below! 

f Number  *  front sideAp _ f Number; 

%  solve  for  y_exitpupil  using  quadratic  formula: 
a  =  1  +  4 *fNumber* fNumber; 
b  =  -4* fNumber* focalRadius_ordl; 
c  *  (focalRadius_ordl . A2)  -  (lensRadius. A2) ; 
discriminant  =  b.A2-4*a.*c; 
if  (discriminant<=0) 

%  if  this  fNumber  defines  an  exit  cone  that  either  misses  the  lens  or  is  tangent/" 

to  it: 

frontsideApertureRadius  -  lensRadius; 

yMarginal _ exitPupil  =  NaN; 

xMarginal _ exitPupil  =  NaN; 

yMarginal _ entrancePupil  -  NaN; 

xMarginal _ entrancePupil  =  NaN; 

else 

temp  =  (-b-sqrt (discriminant) )./ (2*a) ;  %  smallest  root  will  be  the  correct/ 

answer  here. 

yMarginal _ exitPupil  -  temp; 

xMarginal _ exitPupil  -  focalRadius_ordl  -  2* fNumber *yMarginal _ exitPupil; 

%  solve  for  refracted  angle  at  exit  cone: 
exitRayDirection  -  atan (-1/ (2*fNumber) ) ; 

direction_radial_to_exitPoint=atan ( (yMarginal _ exitPupil-lensCenter (2 ) )  ... 
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►  /  (xMarginal _ exitPupil-lensCenter ^ 

CD)) ; 

refractedAngle  -  dir ection_radial_to_exit Point  -  exitRayDirection; 
incidentAngle  -  asin ( (nAir/nLens) *sin (refractedAngle) ) ;  %  still  going k' 

backwards. . . 

interna lRayoirectionl  =  direction_radial_to_exitPoint  -  incidentAngle; 

%  quadratic  formula  solves  for  x_en trance Pupil : 
cotJRD  “  1/ (tan (internalRayDirectionl) ) ; 

yMarginal _ entrancePupil  =  ( (yMarginal^ _ exitPupil* ( (cot_IRDA2) -1) )- {2^ 

*xMarginal _ exitPupil *cot_IRD) )  ,  ,  . 

/(It  cot_IRDrt2); 

xMarginal _ entrancePupil  =  xMarginal _ exitPupil  +  ... 

(cot_XRD* (yMarginal  entrancePupil- 

yMarginal _ exitPupil) ) ; 

frontsideApertureRadius  =  yMarginal _ entrancePupil; 

end  % 

frontAp _ entrancePupil  =  [xMarginal _ entrancePupil ;yMarginal_ _ entrancePupil] ; 

frontAp _ exitPupil  =  [xMarginal _ exitPupil ; yMarginal _ exitPupil] ; 

else 

%  not  using  front-side  Fnumher.  instead  use  default: 

frontsideApertureRadius  =  max (frontsideApertureRadius  _ default, lensRadius) ; 

frontAp entrance  Pupi 1  =  [ -sqr t { ( lensRadius  ) - ( f  rents ideAper tur eRadius  2 ) ) ; ^ 

frontsideApertureRadius] ; 

frontAp _ exitPupil  =  88838B888883B88888888888888888; 

%  gotta  clean  this  up 

%  only  purpose  of  the  ’else1  calcs  is  so  that  you  can  plot  marginal  rays  when 
default  radius  is 

%  given  instead  of  an  F  number 

end 

frontside_Ap_Angle_toCenter  -  2*asin [frontsideApertureRadius/lensRadius) ; 
f rontsideApAngle  =  front sideApAngletoCenter;  controls  how  wide  the  lens  opening*' 

(at  the  front  side)  is 

%  2c2)  aperture  inside  lens: 
if  (useFnumberinputtoCalculatelnsideAperture) 

%  added  on  20101230:  now  we  have  to  backward  calculate  for  the  marginal  rays  from  the*' 
FPA  ! 

%  only  purpose  of  this  whole  thing  is  to  figure  out  where  the  innerAperture  must 
be  for  this  fNumber! ! I 

%  so  we  have  to  assume  axial  rays  for  this  one...  (i.e.  assume  that  image  point*' 
is  on  axis) 

%  fNumber  *s  exit  cone  length  /  exit  cone  diameter  II 1 1 1 
%  NOTE  20110106: 

%  note  that  this  marginalRay  methodology  uses  first  order*' 
(aberration-free)  approx! 

%  assumes  a  well-defined  exit  cone, 

%  with  SphAb,  the  caustic  is  not  a  straight-edged  cone, 

%  but  is  curved  toward  the  focal  point... 

%  so  my  resultant  aperture  calculated  might  be  a  little  too  big,  / 

for  given  fNumber???? 

fNumber  =  insideAp _ fNumber ; 

%  solve  for  y_exitpupil  using  quadratic  formula: 
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a  =  1  +  4*f Number *f Number; 
b  =  -4 *f Number *f oca lRadius_ordl; 
c  =  (fbcalRadius_ordl.A2)  -  (lensRadius.  A2) ; 
discriminant  =  b.A2-4*a.*c; 
if  {discriminant<=G) 

%  if  this  f  Number  defines  an  exit  cone  that  either  misses  the  lens  or  is  tangent ^ 

to  it: 

insideApertureRadlus  =  lensRadius; 

yMarginal _ exitPupil  -  NaN; 

xMarginal _ exitPupil  -  NaN; 

yMarginal _  entrancePupil  =  NaN; 
xMarginal  entrance Pupil  =  NaN; 

else 

temp  =  f-b-sqrt (discriminant) )./ (2*a) ;  smallest  root  will  be  the  correct^ 

answer  here. 

yMarginal _ exitPupil  -  temp; 

xMarginal _ exitPupil  =  f ocalRadius_ordl  -  2*fNumber*yMarginal _ exitPupil; 

%  solve  for  refracted  angle  at  exit  cone: 
exit  RayDi rect ion  =  atan (-1/ (2* f Number) ) ; 

direction_radial_to_exitPoint=atan f { yMarginal _ exitPupil-lensCenter (2) )  *  *  * 

/  ( xMarginal  exitPupil-lensCenter d 

ref ractedAngle  =  directi on_radial_to_exit Point  -  exitRayDirection; 
incident Angle  =  asin [ (nAir/nLens) *sin( ref ractedAngle) } ;  still  going ^ 

backwards , , , 

internalRayDirectionS  -  directi on_radial_to__exit Point  -  incidentAngle; 

yMarginallntercept  =  yMarginal _ exitPupil  +  ... 

( ( lensCenter { 1 ) -xMarginal _ exitPupil ) * { tan  d 

(internalRayDirection2 ) ) )  ; 

insideApertureRadius  -  yMarginal Intercept; 
end  % 

insideAp _ ^entrancePupil  -  [xMarginal  entrancePupil; yMarginal _ entrancePupil] ; 

insideAp _ exitPupil  —  (xMarginal _ exitPupil;yMarginal  exitPupil]; 

else 

%  placeholder,  insert  default  calc  here. 

end 

%  2d)  detector  position: 

if  ( put FPAat Focal Posit ion) 

detectorRadius  -  focal Radius_ordl; 

end 

if  (putFPAatCircleOf LeastConf usion ) 

%  placeholder,  havent  done  this  calculation  yet... 

end 


%  3)  RAY  TRACING: 

%  3a)  rays  from  object  to  lens  entrance:  { object Points-olensEnterPoints) : 
Icreate  lens  points  (points  on  the  lens  through  which  rays  can  pass) 
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lensCenter2Qbject  =  objectPoints-repmat  {lensCenter , 1, numObject Points) ; 
lensCenter20b j ectDistance  -  sqrt [ sum  (lensCenter2Qbj  ect .*2,1)); 
oh j  ectAngieMax  =  asin  ( lensRadius . ZlensCenter20b jectDistance) ; 
lensObject Angle  =  atan2 ( lensCenter20bj ect (2, : ) . /lensCenter2Gbj ect Distance, . . . 

lensCenter20bject (1 , :) ./ lensCenter2Gbj ectDistance) ; 
index  =  zeros [ 1, mimRaysperObject Point * numObject Points ) ; 
lensAngles  «  zeros ( 1 , numRays_perObj ect Foint*numObj ectPoints ) ; 
for  objectlndex  =  1 : numObj ect Points 

angleOffsets  =  linspace ( -object AngleMax (object Index) , object AngleMax (object Index} r  ^ 
nuniRays_perObj  ect  Point) ; 

relativeLensAngle  =  -abs (angleOf fsets) +. . . 

asin  f (lensCenter20bj ectDistance (objectlndex) /lensRadius) *sin^ 

(aba {angleOf fsets) ) } ; 

newLensAngles  =  lensObj ect Angle (objectlndex) -sign {angleOf fsets) , * relativeLensAngle; 
index ( (numRaysperOb j ect Point* { obj ect Index- 1 ) +1 ) : 

(numRays_petObj ec t Point *obj ect Index ) )  =  , 

repmat (objectlndex, 1,  numel ^ 

(newLensAngles) ) ; 

1 ens Angle s{ (numRaysjperQbj ect Point* (objectlndex-l) +1) : ^ 

(mimRays_perObjectPoint*obj ect Index) )  =  ... 

newLensAngles ; 

end 

lensAngles ( lensAnglesCO )  “  2*pi-abs ( lensAngles £lensAngles<G ) ) t 
angles  -  lensAngles; 
nLensEnterPoints  =  numel (angles) ; 

lensEnterPoints  =  repmat (lensCenter, 1, nLensEnterPoints)  +  lensRadius* [cos (angles) ; sin ^ 
(angles) ] ; 

%  ad  hoc  fix  to  block  the  points  that  hit  detector  box  before  lens: 
lensEnterPointsOnWrongSide  =  {lensEnterPoints (1, : )  >  0); 
lensEnterPoints (: , lensEnterPointsOnWrongSide )  =  NaN; 

^create  unit  rays  from  object  points  to  lens  points 
objectlndex  =  index; 

obj  2Lens  -  lensEnterPoints-ob j ectPoints ( : , obj  ectlndex) ; 
unit0bj2Lens  =  obj2Lens . /repmat (sqrt (sum {obj 2 Lens . *2, 1) ) , 2, 1) ; 

%  3b)  entry  refraction: 

lensNormals  =  lensEnterPoints-repmat (lensCenter, 1, nLensEnterPoints) ; 
lensNormals  =  lensNo rma 1 s * /repmat [ sqrt  f sum ( 1 ensNorma Is , *2 , 1 ) ) , 2 , 1 ) ; 
incidentAngle  -  acos ( sum ( lensNormals - *-unitObj 2Lens,  1 )) ; 

angleSign  =  sign (-unit0bj2Lens (1, : )  . *lert3Normals (2,  : )  -  -unitObj 2Lens (2,  : )  . *lensNormals ^ 

U,  :)); 

incidentAngle  =  incidentAngle. * angleSign;  l  need  to  get  the  orientation  right 
refractionAngle  =  asin (nAir/nLens*sin (incidentAngle) ) ; 
fronts ide_validlndex  -  abs (lens Angle s-pi ) <frontside_Ap_Angle/2; 
ref ractionAngle {^f rontsidevalidXndex)  =  NaN; 

%calculate  deflected  rays 

def lectionAngle  =  ref ractionAngle-incidentAngle ; 

Rll  =  cos (def lectionAngle) ; 

R12  =  sin (def lectionAngle) ; 

R2 1  =  “R12 ; 

R22  =  Rll; 
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lensRayDirection  =  [Rll . *unitObj 2Lens (1, : ) +R12 . *unitObj 2Lens (2, : ) ;R21 . *unit0bj2Lens (1,  : ) ^ 
+R22 . *unitOb j 2Lens ( 2 ,  :  )  ]  ; 


%  3c)  rays  inside  lens:  (lensEnterPoints->lensExitPoints) : 

%calculate  hypothetical  lens  exit  points  (if  they  get  past  the  inner  aperture!) 
temp  =  -2* (sum(lensEnterPoints. *lensRayDirection, 1) -. . . 

sum(repmat (lensCenter, 1, nLensEnterPoints) . *lensRayDirection, 1) ) ; 
lensExitPoints  =  lensEnterPoints  +  repmat (temp, 2, 1 ). *lensRayDirection; 
lensEnter_x  =  lensEnterPoints ( 1, :) ; 
lensEnter_y  =  lensEnterPoints (2, :) ; 
lensExit_x  =  lensExitPoints ( 1, :) ; 
lensExit_y  =  lensExitPoints (2, :) ; 

yIntercept_PassingLensCenter  =  ( ( (lensExit_x. *lensEnter_y) - (lensEnter_x. *lensExit_y) ) . / t/ 
(lensExit_x-lensEnter_x) ) ; 

validindex _ ylnterceptCenter  =  ( (abs (yIntercept_PassingLensCenter) ) ^ 

<insideApertureRadius) ; 

lensExitPoints (:, -validindex _ ylnterceptCenter)  =  NaN; 

%  3d)  exit  refraction: 

lensNormals  =  repmat (lensCenter, 1, nLensEnterPoints) -lensExitPoints; 
lensNormals  =  lensNormals . /repmat (sqrt (sum (lensNormals .  A2, 1)  ) ,  2, 1)  ; 
unitExitRays  =  lensEnterPoints-lensExitPoints; 

unitExitRays  =  -unitExitRays . /repmat (sqrt (sum (unitExitRays .  A2, 1)  ) ,  2, 1)  ; 
exitAngle  =  acos (sum (-unitExitRays . *lensNormals, 1) ) ; 

angleSign  =  sign (-unitExitRays (1, :). *lensNormals (2, : )  -  -unitExitRays (2, :). *lensNormals ^ 
(1, :)); 

exitAngle  =  angleSign. *exitAngle; 

refractionAngle  =  asin (nLens/nAir*sin (exitAngle) ) ; 

%calculate  deflected  rays 

deflectionAngle  =  ref ractionAngle-exitAngle; 

Rll  =  cos (deflectionAngle) ; 

R12  =  sin (deflectionAngle) ; 

R21  =  -R12; 

R22  =  Rll; 

exitRayDirection  =  [Rll. *unitExitRays (1, : ) +R12 . *unitExitRays (2, : ) ;R21. ^unitExitRays (1,  : ) ^ 
+R22 . *unitExitRays ( 2 , : ) ] ; 

%  3e)  rays  from  lens  exit  to  detector:  (lensExitPoints->fpaPoints) : 

%calculate  focal  plane  entry  points 
a  =  1; 

b  =  2*sum (exitRayDirection. * (lensExitPoints-repmat (fpaCenter, 1, nLensEnterPoints) ) , 1) ; 
c  =  sum ( (lensExitPoints-repmat (fpaCenter, 1, nLensEnterPoints) ) . A2, 1) -detectorRadiusA2; 
temp  =  (-b+sqrt (b. A2-4*a. *c) ) . / (2*a) ; 

fpaPoints  =  lensExitPoints  +  repmat (temp, 2, 1 ) . *exitRayDirection; 


%  4)  PLOT  AND  SAVE  FIGURES: 

%  4a)  make  name  tag: 
nFloor  =  num2str ( floor (nLens) ) ; 
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nFrac  =  num2str (nLens-f loor (nLens) , * %6 . 4f 1 ) ; 

clipname  =  [' zNO_n' , nFloor, '  p* , nFrac (3: end) , '_ap' , num2str (floor ^ 
(frontside_Ap_Angle*180/pi) ) ,  *_f ' , . . . 

num2str (floor (detectorRadius*10^6) )  ,  .  .  . 

*  _oA  * , num2str (floor (objectAngles* 18 0/pi) ) 


%  4b)  plot  ray  trace: 

%  4bl)  lens,  detector,  and  object: 
figure; 

set (gcf , * ToolBar  * , 1  figure 1 )  ; 

set (gcf, *  OuterPosition ' , outerPosition _ full_f igure ( ) ) 

set (gcf, ' Position 1 , position _ full_f igure ( ) ) 

pause  on 

pause (1) 

hold  on 

axis  equal 

%plot  lens 

dAngle  =  pi/100; 

angle  =  0 : dAngle: 2 *pi; 

nPts  =  numel (angle) ; 

points  =  repmat (lensCenter, l,nPts)  +  lensRadius* [cos (angle) ; sin (angle) ] ; 
%plot (points (1, : ) , points (2, : ) ) 

patch (points (1, : ) , points (2,:), [0.7, 0.7, 0.7], ' EdgeColor ' , 'none' ) 
plot ([0,0], [ -detectorRadius , -insideApertureRadius ] , ' k ' ) 
plot ([0,0], [detectorRadius, insideApertureRadius] , ' k' ) 
angle  =  linspace (f rontside_Ap_Angle/2, pi/2, 300) ; 
nPts  =  numel (angle) ; 

points  =  repmat (lensCenter, 1, nPts)  +  lensRadius* [-cos (angle) ; sin (angle) ] ; 
plot (points (1, : ) , points (2, : ) , '  k'  ) 

points  =  repmat (lensCenter, 1, nPts)  +  lensRadius* [-cos (angle) ; -sin (angle) ] 
plot (points (1, : ) , points (2, : ) , '  k'  ) 

%plot  fpa 

dAngle  =  pi/ 100; 

angle  =  -pi/2 : dAngle :pi/2; 

nPts  =  numel (angle) ; 

points  =  repmat (fpaCenter, 1, nPts) tdetectorRadius* [cos (angle) ; sin (angle) ] ; 
plot (points (1, : ) , points (2, : ) , ' k' ) 

%plot  object  points 

plot (objectPoints (1, : ) , objectPoints (2, : ) , 'bo' ) 

%  4b2)  rays: 

%plot  rays  from  object  points  to  lens 
start  =  objectPoints (:, objectlndex) ; 
stop  =  start  +  obj2Lens; 

%  AD  HOC  FIX  FOR  THE  RAYS  PASSING  THE  FPA  GUARD  LINE  !!!!!: 
x  =  [start (1, :) ;stop(l, :)] ; 
y  =  [start (2, : ) ;stop(2, : ) ] ; 

%choose  maxPlot  evenly  spaced 
nPossible  =  size (x, 2); 
if  nPossible  >  maxPlot 

plotlndex  =  round (linspace (1, nPossible, maxPlot) ) ; 
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else 


plotlndex  =  Impossible; 

end 

plot  (x( : , plotlndex)  ,  y  ( : , plotlndex)  ,  '  r'  ) 

%plot  rays  inside  lens 
start  =  lensEnterPoints; 

center  =  [zeros (l,numRays_perObjectPoint) ;  yIntercept_PassingLensCenter] ; 

stop  =  lensExitPoints; 

xLEFT  =  [start (1, :) ; center (1, :)] ; 

yLEFT  =  [start (2, :) /center (2, :)] ; 

plot (xLEFT ( : , plotlndex) , yLEFT ( : , plotlndex) , ’ r ' ) 

xRIGHT  =  [center ( 1, :); stop ( 1, :)] ; 

yRIGHT  =  [center (2, : ) /stop (2, : ) ] ; 

plot (xRIGHT ( : , plotlndex) , yRIGHT ( : , plotlndex) , ' r 1 ) 

%plot  rays  from  lens  to  fpa 

start  =  lensExitPoints; 

stop  =  fpaPoints; 

x=  [start (1, :) /stop (1, :)] ; 

y  -  [start (2, : ) ;stop(2, : ) ] ; 

plot (x  ( : , plotlndex) ,  y  ( : , plotlndex) , ' r ' ) 

%  parameters  of  plot: 

plotLimit  =  max (lensRadius, detectorRadius) *1 . 5; 

xlim( [-plotLimit, plotLimit] ) 

ylim(  [-plotLimit, plotLimit]  ) 

set (gca, ' XTick' , [ ] ) ; 

set (gca, ' YTick' ,  [ ]  )  ; 

title (' Diagram  of  Optical  System  and  Detector FontSize ', 16, * FontWeight '  ,' bold '  ) 
%xlabel ( 'millimeters ' , ' FontSize ' , 14, ' FontWeight ' , 'bold' ) 

%ylabel ( 'millimeters' , ' FontSize' , 14, ' FontWeight' , 'bold' ) 
saveas (gcf , [clipname, ' zl ' ] , ' fig' ) ; 
saveas (gcf, [clipname, ' zl ' ] , ' jpg ' ) ; 

%saveas (gcf, [clipname, ' zl ' ] ,  ' emf ' )  ; 
saveas (gcf, [clipname, ' zl ' ] , ' eps ' ) ; 

%  4c)  Point  Spread  Function: 

%  4cl)  Compute  PSF: 

%FPA  histogram 
% figure; 

rayPoints  =  real (fpaPoints) -repmat (fpaCenter, 1, nLensEnterPoints) ; 

angle  =  asin (rayPoints (2, :). /detectorRadius ) ; 

binEdges  =  linspace (-pi/2,  pi/2, nPSF_sampleBins ) ; 

binWidth  =  mean (diff (binEdges) ) ; 

count  =  histc (angle, binEdges) ; 

%bar (binEdges, count, ' histc' ) 

%xlabel ([' angle  (rad),  bin  size  =  ', num2str (binWidth) , '  rad']) 

%ylabel ([' count  (  ' , num2str (numObjectPoints) , '  object  points,  ', num2str (sum (count )),' ^ 
rays  ) ' ] ) 

%grid  on 

%  4c2 )  Plot  PSF: 
figure; 

set (gcf, ' ToolBar ' , ' figure ' ) ; 


A-16 


set {get, ' OuterPosition  * , outerPosition _ full_f igure ( ) ) 

set (gef , ' Position* , position _ full_f igure ( ) ) 

pause  on 
pause (1 ) 

binCenters  -  binEdges (1 : end-1)  +  dif f (binEdges) /2; 

countDensity  =  count/binWidth; 

density  =  countDensity/sum (count ) ; 

plot (1000*binCenters, density (1: end-1) ) ; 

set (gca, *  yscale * , *  log' ) 

grid  on 

set (gca, 'xlim' , [-pi/2,  pi/2] ) 

xlabel ( 'Angle  (mrad) ' , ' FontSize' , 14, ' Font Weight ' , 'bold' ) 
ylabel ( 'Magnitude* , ' FontSize' , 14, ' FontWeight ' , ’bold' ) 

title ('Point  Spread  Function,  Optical  Aberrations FontSize ', 16, ' FontWeight bold' 

xMIN  =  1000*binCenters (f ind (density>0, 1, ' first ')) ; 
xMAX  =  1000*binCenters (find (density>0, 1, ' last ')) ; 

%xlim ( [max (-pi/2, 1 . 4 * xMIN  -  0.4*xMAX) ,min(pi/2, 1.4*xMAX  -  0.4*xMIN)]); 

%  use  constant  scaled  display,  always  6.2  mm  wide 

xAVG  =  0.5* (xMIN+xMAX) ; 
xlim([xAVG  -  3.1,xAVG  +  3.1]); 
yMAX  ®  1.5  *  max (density) ; 
ylim( [yMAX/1000, yMAX] ) 
saveas (gef , [clipname, ' z2 ' ] ,  ' fig ' ) ; 
saveas (gef, [clipname, ' z2  '  ] , 'jpg' ) ; 

%saveas (gef, [clipname, ' z2 ' ] , ' emf ' ) ; 
saveas (gef, [clipname, ' z2 ' ] , ' eps ' ) ; 

%  4d)  MTFs : 

%  4dl )  Optical  Aberration  MTF 
%Optical  Modulation  Transfer  Function  (in  cyc/rad) 
fig  =  figure; 

set (gef, ' ToolBar ' , ' figure' ) ; 

set (gef, 'OuterPosition' , outerPosition _ full_figure ( ) ) 

set (gef ,' Position' , position  full_f igure () ) 

pause  on 
pause (1) 

f  =  f ft (density*binWidth) ; 

M  =  abs (f ) ; 

CycPerRad  =  (0 :numel (M) -1) / (pi) ; 
plot (CycPerRad, M, ' . - ' ) 

xlabel ( ' Spatial  Frequency  (cyc/rad) ' , ' Fonts i ze ' , 14, ' FontWeight ' , ' bold' ) 

%xlim( [0, (numel (M) /2) / (pi) ] ) 

%  force  plot  size!!: 
xlim ( [0,420] ) 
hold  on 

ylabel ( 'Normal ized  Modulation  Depth' , ' FontSize' , 14, ' FontWeight ' , 'bold' ) 
grid  on 

%  4d2 )  Detector  MTF  (100%  fill  factor) 

theory  =  sin (CycPerRad. *pixWidth*pi ) ./ (CycPerRad. *pixWidth*pi) ; 

%plot (CycPerRad, abs (theory) , ' k' ) 
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legend {{  Optical  MTF'},  'FontSize ',14 ) 

%text ( 420, 0 . 83, [* Detector  Pixel  Size  “  ' , num2str [pixWidth*1000) , ’  mrad  =  ' 

%num2str (pixWidth*detectorRadius*10rtG) , 1  ^ 

Vnum1  ]  ,  1  FontSize  %  14  ) 

%  xlim([0,lOO] ) 
ylim(  [0,1] > 

%  4d3)  Draw  68  cyc/mm  line  for  reference: 
cycpmm  =  CycPer Rad/ { detect© rRadi us* 1000) ; 
maxFreq  =  68;  ^cycles  per  vm 

maxFreq_in_cycperrad.  =  maxFreq  *  (detectorRadius  *  1000);  %  converted  from  mm*-l  to**' 

radA-l; 

%  draw  a  line  at  maxFreq  on  the  MTF  plot: 

plot ( [ max Fr e q_i n_c y cpe r r ad, max Freq_in_cycper rad] , [0,1], 1  k 1 , 1  LineWidth1 ,1.5) 
text (maxFreq_in_cycperrad+2, 0 . 83, T  63  cyc/mm1 , ' Color1 , r  k1 ) 
text (maxFreq_in_cycperrad+2, 0 , 86, 'Line1, 'Color1, "k1) 

$text (410,0*91, 1 (6.0  M)'f 'Color' , ’k'J 

title ( 'Modulation  Transfer  Function,  Optical  Aberrations1 , 1 FontSize ‘ ^ 

16,  1 FontWeight 1 , ‘bold1 ) 

saveas (gcf, [clipname, 1 z3r ] , 1  fig’ ) ; 

saveas (gcf , [clipname, 1 z3 1 ] , ’jpg'); 

%saveas (gcf, [clipname, 1 z3* ] , 'emf ' ) ; 
saveas  ( gcf ,  [clipname ,  *  z 3  1  ]  ,  1  ep ■>  1  )  ; 

%  4d4}  Replot  Optical  MTF  in  cyc/mm  (measured  on  detector) 
figMM  “*  figure; 
set [gcf, 'ToolBar ' , ' figure 1 ) ; 

set (gcf, 'QuterPosition1 , outerPosition _ fullfigure ( ) ) 

set (gcf , ' Position ' , position^ _ full_figure ( ) ) 

pause  on 
pause (1 ) 

plot [cycpmm, M,  1  .  -  '  ) 

xlabel ( 'Spatial  Frequency  (cyc/mm) ' , ' FontSize ' ,14, ' Font Weight 1 , 'bold' ) 

IxlimUO,  (numel  (M)  /2 )  /  (pi)  ]  ) 

%  force  plot  size!!: 
hold  on 

ylabel ( 'Normalized  Modulation  Depth1 , ' Font 31 ze' ,14, ' Font Weight 1 , 'bold' } 
legend { { 1  Optical  MTF ' } , 1 FontSize 1 , 1 4 ) 
xlim{ [0,68] ) 
ylim{ [0, 1 ] ) 

set (gca, ’XTick' ,0:17:68); 
grid  on 

title { 'Modulation  Transfer  Function,  Optical  Aberrations ' , ' FontSize ' , ^ 

16,*  Font We ight ' , ' bold ' ) 

saveas (gcf, [clipname, ' z4 1 ] ,  T fiq 1 )  ; 

saveas (gcf,  [clipname, ' z4 1 ] , * jpg  T ) ; 

%saveas (gcf, [clipname, 1  z3'  ] , 1 emf 1 ) ; 
saveas (gcf, [clipname, 1 z4 1 ] , ' eps ' ) ; 
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%{ 

%Optical  MTF  in  terms  of  cycles  per  mm  (on  FPA  surface) 
figure; 

plot (cycpmm, M) 
xl im [ [ 0  f maxFreq ] ) 

xlabel ( 'Spatial  Frequency  on  Detector  (cyc/mm) 1 , 1 FontSize1 , 14, 1 FontWeight ' , ’bold* ) 
ylim [ [O, 1)  ) 

if  maxFreq  >  [ [numel [Ml /2}/pi) /( detector Radi us *1000) 

title  £ [ 'WARNING:  PSF  sampled  at  nyquist  =  1 ,  num2str  (( (numel  (M) /2  ) /pi  )  / i/ 
(detectorRadius+lGOO) )  ,  *  .  * 

1  cyc/mm,  increase  ^ 

nF5F_sampleBins  '  J ) 
end 

y labe 1  £ 1  Norma 1 i zed  Modul at i on  Depth ' ,  '  Font  Si Z e ' , 1 4 ,  ' FontWeight '  f '  bold ' ) 

%Confirm  Optical  MTF  by  examining  normalized  modulation  depth  of  object 
fctest  pattern  convolved  with  angular  point  spread  function 
f  igu re; 

f pa Space  =  binCenters+detectorRadius*1000;  %FFA  coordinates  in  mm 
%  test Frequency  =  10;  %cyc  per  rad 

test Frequency  -  1  * [detectorRadius  +  lGOQ ) ;  %cyc  per  rad 
objSignal  =  [sin CbinCenters+testFrequency*2*pi) +1) /2; 

obj  ModulationDepth  =  [max [objSignal) -min [objSignal) ) / [max [objSignal) +min (objSignal ) ) ; 
imgSignal  =  conv (density (1 :end“l) *binWidth, objSignal,  'same'); 

% sample Index  =  round (numel (imgSignal) /3) ; round [2*numel [imgSignal) / 3 ) ; 
samplelndex  -  round(2*numel (imgSignal) /5) : round [3*numel (imgSignal ) / 5 ) ; 

^calculate  image  modulation  depth  avoiding  edge  effects 
sample  =  imgSignal (samplelndex) ; 

imgModulationDepth  =  (max (sample) -min (sample) )/ (max (sample ) +min (sample ) ); 
plot ( fpaSpace, obj  Signal ) 
hold  on 

plot (fpaSpace, imgSignal, 1 r1 ) 

plot (fpaSpace [samplelndex) , sample, 1 k1 , 1 linewidth 1 , 3 ) 

title ( [ 'Object  Frequency  =  * , num2str ( testFrequency/ (detectorRadius* 1000) ) ,  .  .  . 

1  cyc/mm.  Object  Modulation  Depth  =  1 , num2str (obj Modulation Depth) ,  . 
*,  Image  Modulation  Depth  -  1 , num2str (imgModulationDepth) j ) 
legend [{ 'Object  Signal Image  Signal 1 , 'Image  Signal  used  for  calculating  imaged 
modulation  depth')) 
xlabel {'mm  off  axis  on  FPA'} 

%add  coordinates  to  MTF  plot 
figure (fig) 

plot ( test Frequency, imgModulationDepth, pk+', 'markersize 1 , 20, 1 linewidth 2) 
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Illustrations 


Figures 

Figure  1-1.  Two  Usages  of  the  Term  Pixel:  (a)  Pixel  Pattern  as  Used  in  Imaging 
and  (b)  Physical  Pixels  (Bottom)  Forming  an  FPA,  With  Micro  Lenses  (Top) 
and  Color  Filters  (Center) . 1-1 

Figure  1-2.  Cutout  of  a  Camera  Showing  the  Lens  System  and  the  Position  of 

the  FPA . 1-2 

Figure  1-3.  Comparison  of  Ray  Tracings  for  Different  Systems . 1-2 

Figure  4-1.  The  Spatial  Frequency  of  an  Object  Pattern  With  the  Resulting  Image 

Pattern  at  Different  Spatial  Frequencies . 4-2 
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Abbreviations 


2-D 

two-dimensional 

3-D 

three-dimensional 

AMDS 

Advanced  Mine  Detection  System 

DARPA 

Defense  Advanced  Research  Projects  Agency 

EOD 

Explosive  Ordnance  Disposal 

FOV 

field-of-view 

GRIN 

Gradient  Index  of  Refraction 

HARDI 

Hemispheric  Array  Detector  for  Imaging 

IDA 

Institute  for  Defense  Analyses 

LSA 

longitudinal  spherical  aberration 

MTF 

modulation  transfer  function 

MTO 

Microsystems  Technology  Office 

PSF 

point  spread  function 

UAV 

unmanned  aerial  vehicle 
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